“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1988-12 


Bubble detection using a dual frequency 
sound field 


Lineberger, Ernest R. III 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/22932 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


: Calhoun is the Naval Postgraduate School's public access digital repository for 
/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 





























































































































































































































































































































































































































































































































































































ee a ee Se 
iJ ——— ews care rema Pwo + SAMOA 8 ahs Maal Aen ge:aane 
- 5 ataRadpaye 8" ody $e MAME Cede’ uh 8 Meth Uh saad. TpUa ‘el urataerccaretugerannpeatnenios 
s ~ VASE ee », aban ob d. ie t Vo4 6 Aba did My 4 AetcaaRis ceidare 
t ne" SURV # afte tsk LAMA 4 abuees @ BBX pRB 1 AEs, AAsA hima” aOR ROMER Aca 
guns 0r | ; 1, , witee Nn ry fe ree ender aes ante Saad A Ran serteh Alas tds db aise rasecmi@opins rsa 
: : J igh lt ‘ 7 Ft We OSS A, Moat, By lb Ose tea ASA HEC S04 98 saa 
4 4 2 ' aes EM ies WU Tre SRC Ww On eres) : OEM MTU AAT 
' : F +p ; tice Ae oat BH Udo *Ahaargaces Betse Griese. Port dA 
ou ' wege Tae os a ak vate reads ab Pear eh au adate tyqcngs soa Ate,te aon Aang tomrag aaaaes abe 
. Ol - *, oy 2 om Yh PE ron YA a n A.A WAT UM ASNte 25h Peter rte ey h ey en 
° . ; : ' “ ' aE : eal fal A "She te FP te P atte. « fevasd @ Sate ie, . Aike Buvah pee “seri tara a1 AdK® ge, de54 Gn bien 6 1550bAw 
1 ’ Pn am u ° . ee 
‘ a 4 OE " a3 t,? ee foe Shae 0.4 a ire Suits as ate blather dea He ae Gerla ba or a min onaaasennis 
; . , ain? AE Lt ROG Fe be bo: . ATT aM AIM es rhe 42 es £2. Teltuee a’ 
5 : ' i : aR CTy 3 or eg! th Duce dha w yecata a ast: bp plete penne eat te we pier orr vtsasatortgi al iframes: ownpaguapssi enna aa 
; Yass real AY te 8 BR Uff Mdendagialide TET iy, hat rd Pee ccal sasetpa eet Bearaeacagent OU aksnas 8 Nasbesngegnss aise se0e ahs Nee 
i é A t . « 
+4 re wn * . ; 4 4 my i ee a Read \ lates -\asD JAAR ee Teah arte Ach toad ar Rea a ones eee : 
ae: a ee? esp Jae 8G od nt Mehta pat mp heir a ot: wa Raider ara, a few Hit agertsone ie ap a cae SeuHetgB inde BUCgn OP MO BAND: 
1 . 1 ag  . an 4 “aie ya ss oAaged 6 le iti r¥ Wee 461 ASD AABN AS ARO DY ASUS ANA Foc ie erga ihe tee..8 cons a A 
4 ? j 14? jet Mad aN “uh R.ALD Ash. Ar fp Colt 6 eMotie Js, ASARe' WO c MP0 Sag LAR OI Ag: ats 
“ i $ *, : I 1 ees ’ te Myngair'as 0 ad arb My bx & WYER BT AAA OENy eft 4401 90 Frame 
I ‘ 1 1” Rosati alee . Weedegt APS dee mAs ne Ae Aone ag Ww oa rye 
‘ ine Ve 1 " > Mag I Dihe A uhlay Pa htazecd'uc am Ry ede pote ss 
Ned pe Kia) sy , WA ‘7 dl ie mit % no vt pee : see eet Mi iefen @ Baie edieate eee I CN pet 
* t : ye ae F aa bt ‘te ePitser® ett jo Mn BEM a Radner cer MOM sta* wes teenie a1 6 U6 hem 
’ a? t a? 6 @ NMeee gl ene eevmer ye oer) & Aol sd. Cotte Hel aA AS ABAD e aie SVIVIVRT i eh er ee ty} © hgargh ine Bé O88 Lake 
oe ROM abode 1s i TET a aa Reles age ¢ Uke ee Me W309 ott Ge Pages MRS AG, A098 ABs Ap 9 Fb d 144 IAB 0 0I¥ A di woh eh mer bean 
se ee ut why on io AOR Oe oe te ie, as ogng trca® Bak AVAL 8 Masago LS ary rity ak DAN ZCHdeH Bot vagal Ease (Aare 6A Nave! 
‘ oe , F 1 1 a ters " gta ¢ ne .4 ror) be eas Ars niga aw POMS tem ore tars > at» AL ae nse tete 59 DAG ear eh Tee 
. ' es ¥ n 8 adrd phe AA tee Ai : sara Ab 8 Ne “ei Teac m. VIN NZONAG ONS Adee 8 aah 4 nye sean gases 
F 5 Maas e "y wat A he A oe re Brenan iy eS wees Sapa fees one a aime egisin tai tg Aik -pebpbes asm ease 
‘ ’ : : rae yi 1 5 eho - tagor 
- a, ett Be tata ts Fon yt de tee etait ibanez tine ert Mie et see 
' , ‘ x, hee" bt Qotge fs ) Aowrd Here Nia SéGOS 144a8; eh. u, 8 SOARU ats ei 
| Sore | eNO se SN AR Res camer rar ea ec RON 2a ne erage 
’ ‘ 4 . + tae Me 1s: “9 tad ' tgsde “Afndh oe Rages! a Poh tarttidedoa Hs Og ae A Sapa cate gn eabaaas ae erg 
oh ¢ Las Fe Pe MA ede Bangts CATA Oa Ags © Cafe agugas Um alae ge Ady Chm RO Ew 205 OATS 
’ ’ : To. gt b they $0 Ata WMA VOCLAL een 0a Rot ae ager o® A, 7 1 pecans ase Peri) on 
' ‘ Age R GAG Bete 8 AOA the ?; WERiee Ase k ageade Fade 0*R Sele cans bre MAGR9-4 f0.6nger A On ‘ 
' a * tere? Mitt h au“ 44 Sep net oe seb san Lid % Aykroyd nh ‘ ay Or jreae rt oA UEO ag ege 
; \ gt Oona ng teak be jal ain adi ecatentt eit MRR atte endte merge soacaam sen Coste 
i te ' : Z memory ore PedEE, * Vs) BM agticiatty arse wh cect co Gan tek tot? UA Aaees agi gine cw fedaas 9 ee Pp et etd rare 
A ; * 1 y rT ah ot Miss aie AB gM a9 Srmgua sgd.nr beaters tap ngadee sen ao aa; yapaela asad Mies Bee abst ties wane tants 
: a ‘ a ph "tp het y Sahsde igh to" "afr “edaaghen' te ygmeas were aee set oo vane Toliee eras 3s ngcouatne congaane a keine ar we gak AWN MAR pas 
' ' Late pi Pay oe . gah Pe wat Toate aia eenpe Agape ne qs se ‘ : HS sa gsat goto hun bs wed ger yess 
‘ 7, F) mI aay eat amare oeca © A tg a A ad Ost vy : SAjd ene bende URN 
o i | aes oe oa tera! yy * elie apy tpt Bee ed16. 49 Ans Np Me 6d Fa eAA OOESES a aged? p36? DYaML Loe 
' 4 : . tania erate ys yi Betta toate brknnsd Ackeae, GMe4s ATO a “ar Pome aeddet aan 
* a ‘ ; mone la neh % bs (ben BZ Ny WA 10,918 FA 158, 0, Udy + aryh PP IEP Ae on 
7 ’ a } Poe " : » » 8 toa aes tigy 2° : gant ms , Faas oe Besinh sie aittaai se uae de Phaee 
. ' pis £ " j ' 5 ’ 4 VeM te tt § 4 ony i. ® 
: Ac QEbsn PA clay Tae A Ws sa don May * f wages det Wry Medees 6 noceaicales dadrntine Om »0 Cathe 
‘ coal sh See , | i 1 raise often as * pion ; $3 deseles Sou dings gi a ee arscrrecy ae eee . 
-_— | | ' { ¢ OE . mn to 20st chs Ros arrpe are ee ete toe ng ae ze pier Myon bPAbe sy RES sh 
. ' Paes Coat 5 ee gn 8 — te bore ing wh Nenconaté oh : 8 | Pra $e0) Boesoors ta bates . 
. Uh Pdi 8 feb cae toe | s Vee ae > User Ode® Wiog Leetelt aed @ _ a toe . a Nov tae Wiles FAIT EIA, 
xs t Fi spl aia Vea eine Os tt bh Wed LTP gsasaseay Mtl eg \ Sav rat) oes Pe CE OR I cae <a af REE Pees he aa 
. ‘ t 1 , rs : ug cary waka te + Fre RATS Oe yg ote oe Henn BP row mene madera bs ddent! TN, 6.8 ena teevaga Yes tore Sr veer Yd ote ee eo ” 4 
: . Laeh & 7 A ¥ TNE 9 6 Lite Wied 6b Mott fy Juely da rer tual oa peace aT ’ a rarity 1$ tre 0 fA . 6,6 pF IOSOO WAIN O. Bs preset 
ee 2 : ; ; ude UA th; ek yas i £ OTE cece Lhe ; Gone ard-9) 
; : : ; yoga alia og tele Re Ao as tederen tet g <i Beene a Sid ha aA ly \ bras pepe | pr ber pt wep aba i preare ive reson 
2 . ' ' 5 = *4 i” a Poet gywh 8 tt weet 2a, a taro sega ft oiegs (ie A fngethoe Weatdom mcd au ae opie tpehegase STOR HORS h Shy B Oat ate barat de angna 
F 1: ' ign Cae ee on Oat ee res pee OX she gp : DN tat? a 8 ae Art tadoogadt * Le vary want proton de ghee 9, BER CRATA Deeg 85 ¢ 
: 7 lan i? a , ia ts J sue vind aa Yaikrises , : et mM re PR, Bye tele a ‘ot phage nas Pa Vode Fqes ih Fey Mehr iw ot be pbae gy) fee pay ey 
. ab fa : tea Nu ee] ge are tl hy ve 0 Vs “10 a ach fie vy d hea es ane e Sache td joue ay Ayartegazes “a he ee wecrianna Pann oma home Be 
; Pan ; ’ we ‘ , oe 4, «a ; wy Fr ee a mee ath. } Aatned ahs Vt doar be g's Awe! et PDS eat ese te, by66 ter BRA ey due ae Bares Ay taSed anes Oban 
: eo b ' vag 4 1 NY oe Cy aie gods Saghe Rae 4h Reanitagee on aan ) oor were. aek Pe teense atA CoA piolonde Mas emmneragnes 
G . : - ye? : (os ne 3 VE sac any 7 wats oes ’ 00 ay dyad gn Pripesns Qo Ugh ae teas apc Minnsegn sg 
ark ‘ t, . > ° i ” ‘ Sa asitaiee ae seat OCT Re oy te tex Bapave ha 
' ys as ' ; hes Hho vidoes ’ bates wae pe a ee ae ee a Saon Wome es a4°6) iRigivAe apni ta wnnsen vg tans 
US. ’ Aer . \ fae aN Nats pre = S48. or ssctihs nee 4 ee ERD LO Lede agen ry® era tous Meee sogn yaeen 
fi ee ‘ ' ‘ A’ Y Pinon A Geigsergh * 14 98 sanconagaaran shalt ee ale Nees one acanseas 
. 4 i at ’ 1 ning reat ‘ Soke ten Pritt oa hisdicsg eines ’ e. Natedcasasaun T4ee tacos seal Sccaeea 
cen ; : sp ; 3 tasgetty? tg ai tiers Veatm oan 6 tinertene ares: waticecsenl 
. cx? . i * ee Tore t Ae pape se Yelle SA ale he 898 LAN EL st TAY blah 0) eo hohed radant ede ven 
ates F i ' . BAG on ae 0. tet ie Eats ioe a aE Wt oteed coos ws V8 1096 130 9.3 14 Ane Qaifinetas | 
H Cite » Fi # . 1 ; , ‘ rt hay echt e ee 5 ere Fe Me tg Ab zure saga ps0 sorte ed arn aa S29. 7 amen 
. : 5 ‘ ’ ? ' bee e 419s Comte st os it ate % er ee Swit sepona Soave aa ee T£2KW g Romo gs § 
i: ; Dt ar ' ‘ : ti af v5 FN “beh Ret ot eae Sere May dad en bo pe nop Sine enyetes rs 4 ary 
a ' ' . + heme See SNA anes a apeee yee aie Aes anti asenay 8 Prev “it soe fen 
. sO ie 3 erase dese s “eg ir Sgtceeng Bes milena 58 gecvores Sed tp OE inp tae hr. Usdgptse garale youre : nate mites 
ee : ’ - SS ry) i Se th Sere Wyetel neh i 4 os nepata ms Be Pep VYe disgey Sentara Mot ope oom es pag bys rh 
P . o, . wins a yi ; ‘ F os a aise 5, ables pea mr fe aes 48 Sa en Pee epi enee yt > aes tSe Mica 
1 ' oe Te fp igs Seah 0 ee bY d's eens Ve teN Wh Oe fees ae 1h Moo ae tote sone ocrmitie. Tala ge, Rates Sete bed o end 
Va i : : e% HE Whe day 4“ tfegrpero shtenvef NT PEt me Pee iy Barre» trees spt in pane! 
‘ f 5 : $ co : ' o ty ' + tgaa! wt, 4 4 2 wedge LF eS uy a Rae ae then ai eine "3 Dafa ie wr “ta FARE AER ofa 3 4 Otiows ‘ 
. * ee : | 3% ates 1 ate eet e BSS uy ; 
' " e ute on 7 a hpi 8 de Oo tae eS gs ey el ups feat vy * py wp. ent ake nda de? MERes shia 4% ta %9* ese CL pat Is 
. ° . t on. Lae 7 ie { wit Fe a5 tae pe : a ‘yrsP ees 4 A oir + ‘ne nes meen 3 yiag Was karnaas TeS0 Fee 
. ' ' ‘ ‘ wo toad he age, ye fs he Ae) sat a en $4 > at , op PAsurhge, qe" rca ae ine 
pa a Uy En F “ ' ; Ray teey Sey i a P bid eset! 1 ey EERE es Betecel*, ofa * a vant tte iain sys Apannbe 4 “he 
“ ‘ ‘ 5 oo ee og a a | B Thi) ee te Hest ie wy ak At wae he “ts wail g at ‘¢: if f “‘Invin cer 
ot eee : . 3 ers ; oe me ‘Ai 9 MT Anes ag’: A Mats See ASA» 
. nts in . ! ys, 48k attee ’ 
* : an capo 0 Py eae a vo $a Wr Pind on Sesiaet ws Mats sf nqpi Ss se Bae A ‘ iach sapatagce cant 
. a ce i‘ by as 4 . i ‘ ee ase 4 sgota « M3 Re elo hy ARC ST ICY. * i tee ingod pike - 6h .s Sf aes 
. ' . . 1 H ae A , Ss ro + S35 tpt t 149 { fi gatepo ue oe a aa ii ’ aiae Risser Hat Aan jnndtes rie ee 
J ‘ P eek i ‘ i ase is * WSriye To OR ‘ ae " 2 oy ae 3 pete oy vit teeed pase mews ties oa. 
' ' : ; 3 ose Pyare H "s é 
: : ‘ : Mo 943 ts OF ayer deter + ont 1 Meee “22m rey DAS Nan as , a¢ Peon 
, ' x ote oe 20 Re: OR ae) Sane “. do eg pokAccwy 
° 7" Prey "te i ‘ $ * i, ar Att, ce t Shek * . ‘ 2B pais 2 
. va aye . : : t ; a + ay. i . "¢ TE aa Saree tent r in 
A Py ¢ : o gttg us ae i r Phare | Vig are be hidty: 63 ° aa Be 
. *, . 1 Be %» paced ae e uy f a Hi a Nae ah as ie 
é . \ wee or) i! Fi beg ft é ae fie ’ ra 
ee ar: ; . 4 as Ry Ms are lee alae [ogee -taaael i, if Rat reel algpace * 
. as ry oa o% eae *¢ Ces ; , ae omy = . ’ 
‘ tea ‘ b om ca a e th —\ es M nd Aa sae Fi ian t: . Sheets. 
‘ a4 ee eas Pt Gen ra Corte qu ws tat te asa avy AES: rab en “aa Pe * 
b) ' ‘ie ie es B a 4 le LU : 
oe) ¢ Ae Gh ‘ A aN; 4 2 race Ts ik aide 
. H . ! *4 P * $ a ie rhe 4 5 “5 hai f anne wie + TERS ae at 
t F ; r oat t 8 + ie! 5 Fe Bivinc ar iy ‘ 4: Se syleag® mis at 
' ey 8, , ¢ +a : a wen | AGE 9 eth’ ee aia a ata 
. ' ‘ wl ‘i . ny fae “\ a rt “ye tee ETA Gee Yemee Ses i, 
. * ' > z 5 a gk ee 
cd : ear, + i e* ae he Ou Att pais enh, we 
' 3! 4 ae ? ke if “ovey f a » 4 2 Kes fi” % tiga ft o- ont 
. : ' P }¥ vase F 77 oy ay a ¢ ras ‘ ‘ fe iS nfs ? . Rays i © » + 4 ie ai ere 
‘ G "4 eo. . i i fe . . . - * : ortse ate iM 4 me) fms 
: % 1 et a ; . af” 4 { ar ape ‘ie - duets er f 43 te none weak Seas 
0 2 = é } me | 3 . va 
eon or hae Bap adel, us ay oer SENG Bey 
, * ° ¢ 4 “tgsS 4 470.4 
. ’ a ads a 2 va ° f Ta tgsse x Meee 
Z P - , eS at (hen ti Oe eer aL Ts tye gighilecaycs 
y ed 2 a # r P aeert Jes our. ot AED yo wy ‘7 try 2petgin 2 . 
. F : f ry bs bary ; = “he, op '* lew ee fs os ee 4 47 oe o Ahaie 
: : : : Kes Foie H8A 5, it | . ny th, ar Bee eee We ie 
P ‘ ; ; a ae o sf A ware a st Re“ (i i fines eit Z 
: . ’ t _ Me, etec s of *: i 4 - 
’ : F e ; : cs ' 7 ; . oy ae " “y Ae SE 3 Ree ym ' ee vey é . st : sei! mete! A 
ne . ‘ ‘ add Po U i had Pyas af ‘. ry eens 
' oeuk : 5 oe gee Ay Nan Gig nt. pg igkeng, (es 
= ’ : ' ‘ . = G e ar 1 ah} Men Ae ae / rte oe RR Fi 
. - st e 
. Fi , ues . ° a O46 
’ 
t i ‘ ? = fe eae 
. ae ‘ ; : : a Pi 
po oae8 ’ Foo 
. up . ' a ’ e ri "eee Nits Seth 
' wie $gksn 
, ' ‘ * : . hate fy: 
Oe epee 1 F &. ve ite « Pure 4 
' . Oh a ' ea a i} ’ 
. fe #. oe ee oko ae Pay 
. ’ . 7. s bat ” - 
. ‘ 4 6 oe v aa, 
ne a L “ ‘ ist : 
2 . a ' 
. 
. y t u ~ ° %, z . 
‘ . rs = ge nce ‘ : : t at 
; fecmoul pie 
’ ‘ ' + owe | > % 
. “hoes ' 
Piet} De ’ EAsre rs At ve see fF = ni so 4 ot 9 ce ores 
* wt . § = = er R mets Cd ’ 
j her aris" ercrsst th é : Se gts ya i tame — ee A Yeeeah eet 
= ° F F ‘ ees fhe nT sald 4 res Pee ae on yieia iit > ie (EE ES a 
- ' hays ; rut yee nee fi id A <_ $93 Kyi 7 Ey a tye rs CGP eyties 
i ' : is . id 5 MOMS oe has 7S Cp, Heer: 0 As 3% eos Mes “ 3 maby other 08 a 
3 eK phen a e a Fo] 2 a a Wt 
7 2 ; Set a. a bei p ad ok Lean fae Le gees aheta 3 iets rites trig 
Neen Cae ai ‘ fc rts bea ce OY rl ot dss Us Gh Aqine (oaatet verdes ined ride Eo tron eg a 
: 2071 wae st poy ee Ra ey i ae a Sts ad 
aoe , : arnt ' ‘ aad Co pamees eee # ue os 2 ber te oe ae ob ep ie seek he ie 3 
a ve ‘ og 1 ‘ : He ptt i pe fi Yew a here one r ay * ; Rintey Ls Bae’ reas 
° . pate aoe . ae x as eee gat: Fe ts eee ~ Se “o-. You Goi Wen we se ty . 
: ‘ ‘ at ga of tetare oe " Cr al Re HN . xt 12 Af Pe cm +aghe, os syhtete sou: 
. ‘ ‘4d ’ 2 nines ge AMS ee ang VEN, oth Aad or ge ee oN ak POM 2 Maoweie® 
P rs A r pe oF of ree , ae ey ine: 4 x “9 ye oa 
. ' : — an Aa wet ot vite ; . 
i ‘ re a “t * . 4 erg # fof sve MHelfitety quit tc i ie eee tA 
“ ‘ *. 7 pie eee obey ot eey On 5 tt Kas Oo ah # ae 
; t ‘ gee z P ‘ a *s - Ri herds a ah haa pA coe i te 
‘ ee me . yo ca, tanle 
D . 4 j , ie ig : ah A zoehe , ured Re > yy Hes 
. ul . 1 ’ ¢ * ' ; F - , pees ray of i Ot. i i z bya ig geen 
vo . > . : A . . ¥ jt a . nave ene . pee tee 
' . ; wie a . 7 3 F ‘ Ras Pade a oe " “t ‘a Gen : Pr Aen Se : Gees 
‘ U 1 s , r yw afha . 
Fe Dee RCO an! rE ‘, t feds dae’ ae ak LOO PEE eres tes “Rane vA, J ee eae 
. ' ' F °} vi ' sty ety of oe A ime v ‘ by z aig add 
: ‘ ’ : 7 pes ‘ke coy eR ge . a6 hai 4 Jae fe ae fi a a hy dna, re ig, ; ¢ a * «ee Sota ip fo nin As Aad een are e Bebe pane pres ret 
A ’ A . 4 fe ; nr ' © » ae t at wh pp ieee £ A Hy he yous 
' ' at fe ry ie elet ‘ ct é ' fe Ais ae is oe bas a ye] ° e Ay ag Bis 0 
Fy te ’ Bl F » t ¥ Pret) ( ‘ ele: heh = bp a4 4: + e pel wih 
¢ Ke “ sale - a + * =: ra aa owe Hi os 4 fi, ; BP ¢: x a. tes negt mvs wn q ¥ ay HY Agta ie 2 w hele eagle: cerwttery 
- . . Ge rare. one : mage ‘ ; Pr 5 Fa ota ies Sebi ook ye yee ieee af AP: ver sPevran 4s 5 spb ate arias Bene tgs 
Ye ain j tees ghee ey Oy a ee ona nc oi wd sibiped oe “ae ie ior amet Seestgty pdapbasey gabe muni sing ne Soeet 
‘ “s hee y , ots + Oehy Are Gar he ee ais: e: a ar ads 4a gus 6 YS ch oli cae ath ene © Spe Vigyrecua gies pane Fa ie Sl oareat dd 
; U ' 7 ove a o3 ‘a ‘ Wey 3 ue ke Ayee 4.4? aie r 17 ar sp ee gt ae ; uy bie oti, aE, 858 pigs ta emer tale ye ay do. ep sae 
’ Me ot . mans ts specs ye oe platy 2 ie ee x A ba wala verge! 0% 4 i fas of ot a pr aretate Pr MTer CTT ob! ‘2 
Aa PP ‘ ’ Oye a of oh ' g 2 tte @ ai YC Bae nite red Bera) seve tance Foe tam sp arypeoraves 
2 . ° 1 ie . tk ' ‘ at ” oe? asf Yi . Pe why ine Fa Seah Abe 5 é apc eee 73 RE fa ein re Atom oes i ie 
Ca ' . . « ‘ . Pate te, aan ar ue uM _ a psy <0 GMO 20 Leg gh gyn ge 
ee i Seer we ee iat SP CE aCe LET De EO SY MOA cw srt pita ented eaenagen 
te 5 4 te Pen er ee ee koges + ‘A es fee? Se adel Meee ob ¥ c a, aye AP wm on ape 
' ; ‘, TO . Ub Pg Deal? er ‘ Le ° - se steutamer o 7 1 987 pemeuih aie oat a sarin x " ar eb ph ART a gh ahve Ee areas Weaemce acer Pio de seagone oe mafpiety Chsttatigoretoess - 
ee ’ ® 4 af) Put ae 5 Be ©, eg r * a > Boe pore oO of 102 8 p42 FYR OSD 6 ® eee ae ea ‘ tea! ‘6 Ge yr Gary G waned 
; . oes Tee Gee Ey fe eran aa Oe pmte oe tae onion: ae a 5b 10 ,4e8 ee ie hy ia Tin" 96 call tA "Bihan aes sree ele levyarte rs. nid cvecniet eas par agp ai 
' . . o4e ‘ ‘ ‘ 24 , Ws & Uh Use ek ho a , a@yre a , | Sew 6 ee easy slut ge, 8 on Fhe 1p veut bee Seb tere ame ol "# oe veast huyhemeae tear 
en ’ , ; wre be ay P Wo Z ag = ee ae OREN ae eA oh p¥e Ne pe ah the pe) Ue ed oe ety ce wigs sigh ee en na hea ety Ores 
, ; ’ Aa eG ina Par a ae ob fre eed , e me we 4 & Lape « ore vA ese s <. Of Tans § lee Ape 8 FUND. ye nog Ce he goed DT elder ee riba jsp Yn) mated 
' F ears ae» LS oe: Fe ey BO ie ee cute) OLY Pech $08 ' the eg es ee via! hye pesado Piel pruttgrae 4 ete zs Z nangge ornereene teal ynts on mae 
: ’ one oe os 1 Has Po i : i Sti pat she a hae Oates # rear mee PLE | set tpl A pte er apes tom ° 
’ 2 vie ' x . “ ry iy, + 4% tee ee Sas we Bin ota ee . babi bhak pasar Akaaie eee he 
mer t ra ra nan tree! He s oe ats : ate Ais = a ey eee as v Dis met ed Fi reales Pane iene ' ps gol hare tpe ol co pubes Cefn erse lc Corpse ane Lt, ents 
re arene 4 a a> J cts foe t@ ate 7 uf" as a7 eae ae erate or : “ses fir. wi sto? 2 : Kas aiueted> saletiboe s AFF pe TAT ees . a“ Pps . 
> 1o8 } aad r . rf ae a Fy ro PF we elon ae of BAB i be | death 
ae : cone Ast eee ae SPC iar. MMO Sete. ys Be id Bete ¢ EE ee se ba eal } Pp ietngrentritne een Seat cerns 
Ce ¢ — . tgs i » Ld 
/ 7 ly en Ee ee one oe te ee y 2 FP p itn fatyy ot fi ariosyi Sete Sle ass atten ee Bg at ntae thot cater reveyte Haale pe sy nsntr 
. ' . te L P ‘ a fi ania Os “4 ie . your ae 
: ¢ So as os ‘ , ane ' vies ie vere ‘ 4 5 4 ne nes rar ye ‘ 4 ie $ “ n> patients ee a Dor ‘neo ‘s ad itetinericae apr Eas fla raat cath eee pe ye elys iyrsstuneqares ocoe 
Ug . . * . te a4 ~ 4 bam ge tyes i ety * Ui i aa 4 Bybee. Lato" ae Sounds & we A ote OF wd rate? ehia, x . as orga 
6 ' 4 ° ‘ wie fs 4 s yea Bag oe tha ty — wy af a Pee . "L aSate ee Wm eelinryae gt 8h Wdriegyia OMEN, oxTee thd Col sh OP eek hated ad aed 
ran Wiig atk ae aoe re eee ’ Fi ” wen et eee a 7 i ete Be os eebicQoiees : 14a 5 Saele! 4m yehas ilar) abt's Pere ge Geen’ " a Saea WAR Gooey 18 TKO 
' " " US ae a peer ees f fs gee GMA de Spey Pe rf fa pe *y tg Petpet oath nye ae ees G Taveds 39% otys sr arate ARE o ANE eras 18 seyret stem 
: q ‘ DOs if Bees ease iy el easel ous “gs eed Fe tee A eek aye beby AECuABIENS 500 ne Wane A vl Gabe 18 yer th 
‘ 1 7 HONS Uys “Dotepeti ec ye rhe tbe ng aves tots Ny mips rer ee Pees wh Sesenek: yw ie Me are macaw \ kn Petree aco faery rill area 
Ae Lye - ine : siete ot ye er laiencs lla a tte ri Ada net le pve nye 1 A hth 488 Sha “7 *s PAN ‘awh eb mouiinte Plt Ab mal ghd Speman Ber sit Peptide 
. ’ ‘ Po oF & 1 5 * =e a ee 2 vip 4S cele IO Arde why? oha a eden 8 Pyare PH ges e. git) on Dewnty: graye had vty ; oA: erro) 
. et ; of . . . WM m- octets Ff Se a Lt bef ole . . Sires a os i 7 “ geyper Pept tate ge = ” ‘i L oe ah tial sf wt etyesigten. 
A ‘ Cena i te > Ses i at} 5 ; Fe shore . . Fate yi onerens ie ar (a) eal aetna eyed ewe wast eve 
open : ' ’ Le id abe > 1 i" ti ° 9 nti eam v ops pres harvey we " belies he ri Pee ae a) S 
4 ' ’ . (2% ae S i) are ape et . ete pia ne ke rene oF: Ti RP FAH 
ae “ae ‘ Z Be aA So G ea : Ae saan a: caer = oe a Sophy Mu fie ae ter Tiss . hee 9 wren ey os Pcaby an abenae arian ati ane an EP Yet st rivsahry all g Se 
Jee ae ke » 4 Pes iene : Tg tenes ig 3 os mays eRe i) Sapiens Wt: ep peta ster Pb te aeiats beaut te Sens paee aoe veg ews ienen pp seat abe me 
1 oe " - ae 7 , io ~ +e 1 Tew mpeg tg tee \ its re airs a 
, . eh te a Wee J Lo Se | ~ Lied DMO ope bIe' key a a | 9 va Tupene tat a" ri 
me) SRM aE Ye a Brie RL MO an br eae ag aero iden Aven ar et a eynne tate ware 
. F re } Eniege * : saa. oe ry Paul Odo aM 7 rig dysp lat | te ar es of be: {aed white poate see B a Dvert ee iit apne an id ae 
. z ' x > S aa « ‘4 eien 8 = pee ses ‘ hp WL 4 del ss: WPF oa = TpAgeS py genenseswe +, egcadrdesds et oo yt s ara Tebat Aver an is Wo 
‘ y ” oo! LE sl al oar o 019? : eu ‘ Te a a eae Pye a's x ’ 2 So 0a ats Aateng 1 a te weg POAT ONT ZS oe Patt brepilaig ead FO pToer ow rQne 
1 F ee cle rs va aie nes as " Ayre os is 4st « Bey HF rz eae A yeas ane mee aphecte e 288 aly Vo ravenve ase ak a eo 
; . +48 ee pea'y a eee : Dari ad Me tye) a8 babes , Hae fer, Paid) sitne ueerater ot Aon ame: ewig rea a The gatet yr rw PTTIN TEM ANAL a Oy Gt 
2 _ 7 aL ve 4, A a cre te ae hed oe gp fe  L een Ye Ribti. Gn a ae i i yeseen rae § panera ae Pets. reece eset Hopteweine arin 
, ‘ an Po ey visitors et OME OAH, WAP 9 ¥ O86 pO eran Feta? see! ‘ Doar Tl oad ie * 
' = Mune nee "an cana". sim ONS my ibe serene _ ote 7 AN a pres egta tiny A 5 oi Zg0 outa fey 2° Sige 242 OF a arate 
re wr) ' ae 1 ' % OG a Rin ae en a fhe. oem i av ot 4 . of) Weale® = 5 SEPP. WAL i i eb gs 3 Pe ae eer 4 Feet A ahha eran esa a3 og Pek es {OP Ney oae Spree won 5 gon 
i ake P ” 4 ; wae -ae - 0s ey ie : a fe A sedans Pay Pee. § a bie WET Dar ee ik bp ie bk as phe? ties Tudhyedsa ons eee EES! ere ee Letty ies 
. ‘ eas ' me id es a F< ape 4 ee ty AI Tr aD ER 21 i exengratanhe Ae 2 gta Te bay Joga Wp galeale? o8'a ge 5° pF if FOR Oy é 
. é os mre mre oP i ace yeys CRAP es Pers on ages ooe phn are Toe ky sh rete La oy Pe MAINS EM Ye TORRE? 
in MEE ae BS < Vorecesitn ae ee ; aes F tower ‘ Jl Peer ak mr, Hea See ae aati A rs re “% tee? Boe Bs ofa iil sete ns atier! titans Mr eee: kere os Aqsa SM hohe LB fe 
PY o of oh "ye e Z . . wigre os wie 
’ ; 2 ' Cae AS, rc oe a) tee i baie wi ie reid og vey Ate =e fy ost Seema det nit gee AMA aur: ag ky Nelioa sl gash oi Pe 9 Fei B 
' no é t . Bhs te ‘ - x Jee wht es anny? Fs : iy ate su yrs v° "wre of Ts er ae eens 
! Tetoks cme ; ae ve ane Ls) ee : ce * abst is) eapsgt age. 1” x bee weasel HN oni ates fegetontcas ie oie eis arate ynaror eter cet 
% F r 1 1 ' tow ‘ vey 4 “ : a : Hi 3 an spate wat Ol? Ate CUS 0g G8 oom yeh osreb MF IF PLT othe 1 ttre oe ts wae eet Mee ete 
; y . 4 fd : © 93° 7 ye i i a Shee 2 VOR Wa, 
Gn itl F ae fa in sf ph $4 ob fe bc ae bgt wrt he A area av : 7 oertary were prota? SGM alts There VAST My ory gt a ROTA ok He 
i ans a, eee ats ae % 1 “3 ' eat Oia tee Gea Oued . seopheyt Sif aren tee OY oe erie cai 
ce : re Se aye ‘oe ti x ae "Rt a , ie rt (50K soe > ee yee Ms oe wae An tere garter te 0% at stay Ride ean peed pe ae Mier wamansints ot genera 
4 ft Pater one a tay ‘ ee 1 Pe a #8 of be Gat sete y ee tas a, 6 fou HOG ‘ety i siaatasaenes v0 tore the 16 wr Polen ymca tbehet Maraseegaes eo eebotkrgts on fies or taye 
te elie ee ‘ , : Ya 2 estan ie ovEGt fer 8 TS ed aes okt a pe od Beets a asrelnts > 7p 0 inate REEMA TYIOO «1tm Ze PN AEA Te oor tasat eRe 
i ; y j 1 84 bate a tee & Mey aig ayever Fn 08g ee 6 OP TL ek, otk Go cle da Vin Lk adie’ ob Bleek ee rs Pd uy 
0) aie ‘ =e) ' s “4elos ° 6 él ! weg "I aie 9 ool saat 1 (44? An etoe Ve ee ee Mey OT hb Ra ae a ora tLe ie 0 AY. OP relent aie 
, , . wens ate m4 4 .. os a eee Le YP e aPmLM ART. Ca be Bt ig Moe We tea oF seth y Sac an Nd hy Sy he Re URE aE tae OOP uh pevsers MArphine oe oho 
: e ss ji aa a ‘ ae nw re debe 52% See uae % yb tek 74% ] Page Feast eure s G87 oohitemet ot ele g He 4 Svyhon a Syae Wr vee va * ee bed cise: 
' ee sae ¥ Be a ‘ ‘ a ie a ves ae PO it ; a& _ oor mnt un wate Oo ee matey age SN tata aie eens cosy Gevevtote "aes, 
= ‘ ’ = U oe ge es eee t ea et bretnte Lamp. He Putts me gees nym eH teqae : ar 45004 q qotm Ld fire eas 
8 ' (Ue ae ’ ts L os FT eer cme Pa war 2 od%ey) 2 earnest ms Be or GN a ee pegtres: eresty' ., CW) be ht aaah et) bi. that a poet sreeatourense 
oe a S sie W ie re . f yn) ppl oh *. : 74s he 58 weary oo . Cree y. % gepared py 'e'# 4) eye, Waetee hm Spe ve eee eer Ores et ge dighaewwar 
, on : cs pe Se Pe bar pant Sas ab FR eka tee Sevier ey ota wheter Pela Pigtabae ge Senn CON ISLit ee TROIaTR he a 
“ ' _ a Ae A UE hal Oy o, ze e : ° area » Wier ger pl oh eld ee roth we 
é Seerece? § r te ee ey ee . ares M1 ot at deters h a et Lacouette ner emesert - vt eddie mae ae prac ere erHe ITT. og eH 8: 
1 1 ar) : 44 8 0 ne faa tebe a. cor went ote —aledibiass 4 pple Me lirl . she/arrhe Liab eye hae ce 
. f ‘ U + ‘ r hone 2 cM be ae 18 extern 30h Tabet) io EM nda Dp Ye ih patna epee GO" eg unTmT Li Sapna: Aenean 
: i 2 ‘ ai a “ an Lie eek BEd hed 5 Fk ghseskpe one oe eer i Nae St IL he dpdpd aa ag. 
F ‘ ’ F ‘ ‘ ’ ' oh WN we Naleer Cas UC) es he a ft my ght: A tebe tan ek ee NOTRE OTS 4 wr ARTO "88 OF wee ty 
. . % . * * 4 een ew = SP ote om hare eye ‘ 
. wre eel iee ‘ ; A a0 : wey aS ns whee ere es tta st ie usapcay 3 . : Fviig ay leealyae eieteunne + Este eed Bites oS CURRYS ao ew ; anh capa VAhIN VP omacbe ee ences 
‘ : fe I, ‘ t I Baa he e Erol sips i TI La Pht : Sy" REY 4 Me snes etc fag: ee a Reprrtt Sae ear oA CeO PON TD OSD ew ys 
CAs as pa 2. eu a! a wipe 19° th Eat! 7 Berk oe mi ae iek eat oe aber! ‘ 
, “en Mg eo el EB te bet iat ae age Colom em es at cm Ryding oa Bere Stren awen! mere ppeeieys ‘Gee ee 
patter . U s ‘ “ ‘ ‘ i « WETIA a5 oe a @ A 4 SL a ny + he 2 pS niginl wre fe PePewigvertey Igemeurer es ets as rere 
‘ ’ ' He a le a 4s phi eae xs, (ase duty’ : ne hy Sirs et aca ware mati. sis (at egnenaes Sepia lhe pe rere ine wen 
. i F ‘4 7 ay on reat uv ae re} H a, peas aE 2 He an cage ’ i] 34 MPAs ye CO cay path ie tek cok be B 
. ‘ wh ees o' ; te &, 49 | eae i ee Fig v ee oe , £5 wncuinbyerdeeiey Vee EIS ow GMP ry thie SEF M MT Ge gs greet = 
ou ¢. i ‘ Ce b ‘ é ve ae r ig ae | i PA ay {eres ee oe a rly 4 rer 1 . ° ination vin -¥ aaa ne A lye awe g rp ig Awe Oe, 
oe ‘ Py Ci ' o me nee iD poe “die td & LBs PI CTA Ue py Pitta he sian 
“ 1 oi r) a ie) Ce é ot. Cetgreht pial = ba wn Leda td 0 oe a: ae pia virer rains sues, 
ae ’ f * VALS Kore ¢ ee . 4 + E's sia et fa dat faa. gy oeeetse 20 agp nboLet GPF LS, le VA Wr AGTA hh a 
{ ‘ ; ’ cy i ‘ity ” 2 ‘ats 97 kt" Meron Wehig* age ment ey ree At ° B, 
Psa Uy my Ue ri aR %, bd torventing TREE oruedige nate es ee at a Ra a ee pear pacs Gy 
i ‘ ‘ rer rk + - a tvee as KINASE or as is wy a wrens *: 079 + ITD ATH NP beet, 
. “at ; 1 8G 3 i ore iy : ONS vs rb o Mess iy nee righ ms nage LS b atghe 3 ® lh uA) fe ik Fs a od ‘ae 
z a <a is re e, Dial ities l 2D te loo ee ¢ Fa ee TT ey Ce ee te ate agi yehene eeaie 
2% % : £2 API PUM OY WE Guetta ttn go . a yretnbretel ee: ‘ PTT Lis Sad bee 
é . ’ a : A : 4 5 are A ay ¥ tiilm ne “ fom ah ews Teen eet tel oa? BAe th ae Rea EA A Ae ean 
ae, Du = uy ns ty uaer TPS oe Par ros ¢ iq Ren el : sees arate h Fa 4 no Dade awh mak PER 26 AH oe EVOL Oe He OTs Hy Seen mire aU iBT G 
. , ey cr 8 gt sor et * sa" ise factve! MP ye Spgeypte eel 8 ee” aah 
* . ry . 1 ' “ty ‘1 * r roe $74 . ‘ + “ate. 4 i * ; MT SY aL Gy bal 2 
4 . * . Orgy yi 8 ' = e 
ae : ‘ ' ‘ e Lr eS -_ oh . i 
aes 
oe 











NAVAL POSTGRADUATE SCHOOL 


Monterey, California 








THESIS 


BUBBLE DETECTION USING A DUAL 
FREQUENCY SOUND FIELD 


by 


Bunest: RR seainebergenr, sti t 


4e¢@ ¢ 


December 1988 


Thesis Advisor: Anthony A. Atchley 





Approved for Public Release; Distribution is Unlimited 


7 ese See 





lassiiiea 


Y CLASSIFICATION OF THIS PAGE 
Neen hte a ht el 


REPORT DOCUMENTATION PAGE 


YRT SECURITY CLASSIFICATION Ib RESTRICTIVE MARKINGS 
lassified 
JRITY CLASSIFICATION AUTHORITY 












3 DISTRIBUTION /AVAILABILITY OF REPORT 
Approved for public release, 


Gistribution is unlimited 









LASSIFICATION / DOWNGRADING SCHEDULE 


YRMING ORGANIZATION REPORT NUMBER(S) 5 MONITORING ORGANIZATION REPORT NUMBER(S) 

















iE OF PERFORMING ORGANIZATION OFFICE SYMBOL 
a1 Postgraduate (if applicable) 


ool 61 
RESS (City, State, and ZIP Code) 
terey, CA 93943-5000 


7a. NAME_OF_MO 


ITORING ORGA 
aval Post 


Braduate school 





7b ADDRESS (City, State, and ZIP Code) 
Monterey, CA 93943-5000 







9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBcR 


iy 


1E OF FUNDING / SPONSORING 
ANIZATION 


8b OFFICE SYMBOL 
(if applicable) 





RESS (City, State, and ZIP Code) 10. SOURCE OF FUNDING NUMBERS 


PROGRAM PROJEG@I TASK WORK UNIT 
ELEMENT NO NO NO ACCESSION NO. 
E (include Security Classification) 


BLE DETECTION USING A DUAL FREQUENCY SOUND FIELD 





ONAL AUTHOR(S) 
eberger, Ernest R. 


DE OF REPORT 13b TIME COVERED 14 DATE OF REPORT (Year, Month, Day) 415 PAGE COUNT 

ter's Thesis FROM TO December 1988 128 

LEMENTARY NOTATION The views expressed in this thesis are those of the author and 
not reflect the official policy or position of the Department of Defense or the 


ted States -Government. 
COSATI CODES 


_D GROUP SUB-GROUP 


TRACT (Continue on reverse if necessary and identify by block number) The design, testing and analysis of a 
al frequency system to detect and determine the resonance frequency (and hence 
e size) of bubbles ranging from 100 to 7 um radius is reported. The resonance 
equencies were compared to estimates based on the rise time of the bubbles. In 
neral these comparisons agreed to within five percent. Although the system is 
t ideal for field measurements, it identifies important requirements concerning 
e size of the sample volume, the frequency and amplitude ranges of the sound 
elds and signal processing techniques to make an effective system. Bubbles are 
stinguished from non-gaseous particles by the nature of their nonlinear response 
the dual sound field. The system is versatile and can be modified to suit 
ny research purposes. 





18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 
Bubble Resonance Frequency, Nonlinear Bubble Dynamics, 
Bubble Detection, Bubble Size, Dual Frequency 












TRIBUTION / AVAILABILITY OF ABSTRACT 21 ABSTRACT SECURITY CLASSIFICATION 
NCLASSIFIED/UNLIMITED () SAME AS RPT C] pric users Unclassifie 


\ME OF RESPONSIBLE INDIVIDUAL 22b TELEPHONE (include Area Code) | 22c OFFICE SYMBOL 
essor Anthony A, Atchle Code 6lAy | (4028) €)46.22)8 C1A | 
All other editions are obsolete ay ab ee ao 


i US. Government Printing Offtce 1986—606-24 9 

















(i, 


0 
= 


1980 


i 
Navy 


ie 
S 
LEO maar 


~ =i 
Saf 


dad State 


of Southern 


by 
sneberger, 


i 
Ors. t 


Re 


BRnest 
LIGeUtenanc, 


Universi. 


(L 


Pal 


~ 
= 


Ucmiited aa: 


Ss 


ACOU 


ENGINEEFING 


~~ 


Brett 


- 


ae) 


= 


- 
v 


Gy 


wAT 
Noes 


pa 


nn 
ee 


en 
ey 


MAS 


Fed tind 
-y 
“ 


ete cx 
Ree Ie 


6 Ve 
a’ 


ee 


ta steed 
_~ 
1 


Soe 


% 


=), oo me - 5 
So. Lol eee 
ee ee od ~~ hed rs 


a 
Sr 


fr, 
fr} 
E-4 
Soy 


a4 
~ 





ia? 
Oo) 
Vbe 


oa 


of 
Dp 
te 
ie 
DY) 


ths 
ba 


—————_—_—_—_—___=_—_—_—_—_—_=—=—=—=—— —_ = 


ABSTRACT 


The design, testing and analysis of a dual frequency 
system to detect and determine the resonance frequency (and 
hence the size) of bubbles ranging from 100 to 7 um radius is 
reported. The resonance frequencies were compared to 
estimates based on the rise time of the bubbles. In general 
these comparisons agreed to within five percent. Although the 
system is not ideal for field measurements, it identifies 
important requirements concerning the size of the sample 
volume, the frequency and amplitude ranges of the sound fields 
and Signal processing techinques to make an efective system. 
Bubbles are distinguished from non-gaseous particles by the 
nature of their nonlinear response to the dual sound field. 
The system is versatile and can be modified to suit many 


research purposes. 
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I. INTRODUCTION 


Bubbles in the ocean create a wide variety of effects on 
sound propagation including absorption and scattering of 
acoustic energy. In order to develop an accurate model of 
these effects, it is important to determine the number and 
size of bubbles in the regions of interest. This thesis 
reports on an implementation of the dual frequency method to 
detect bubbles and determine their resonance frequency. | 

This implementation was conceived with the following goals 
im mind: 

* Detect bubbles with resonance frequencies from 30 kHz to 
400 kHz using the dual frequency method published by 
<i (Ref. 1] and investigated by Hampton 


* Measure the resonance frequency of bubbles as they are 
detected. 


* Distinguish between bubbles and non-gaseous particles. 


* Determine the parameters required for a high probability 
of detection for bubbles in the sample volume. 


The dual frequency method is explained in the theory 
section. Supporting theories regarding bubble resonance, 
visual determination of bubble size and nonlinear mixing 
properties of solids are included. The problem approach 
section details the experimental apparatus and tests. The 
results of the measurements are then presented and discussed. 
This thesis concludes with a summary and recommendations for 


future work in this area. 
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II. THEORY 


The dual frequency method of bubble detection and sizing 
relies on the nonlinear properties of bubbles undergoing large 
amplitude oscillations. Theories concerning the 
characteristics of bubble resonance and nonlinear signal 
generation by bubbles in a dual frequency sound field are 
presented to clarify the dual frequency method. Also, because 
one of the goals of this project is to distinguish between 
bubbles and solid particles, it is necessary to understand how 
sound is scattered by solid particles exposed to a dual 
frequency sound field. A discussion of radiation pressure, 
streaming and rectified diffusion is included since these 
mechanisms limit the allowable pressure amplitudes. A section 
on bubble rise time sizing covers the method used to verify 
the resonance frequency of certain bubbles by an alternate 


means. 


A. BUBBLE RESONANCE 

The fundamental response of a bubble to a perturbation is 
a purely radial uniform volume pulsation in which the bubble 
oscillates about its equilibrium radius while maintaining its 
spherical shape. This approximates a perfect monopole source. 
Other responses, harmonics, may also” occur but the 
fundamental is the most efficient radiator of sound and 


dominates our interest. 


The frequency of this fundamental volume pulsation is 
determined by the stiffness of the gas plus the surface 
tension and a mass contribution from entrained fluid. Many 
derivations have been published relating a  bubble's 
equilibrium radius to its resonance frequency. One of the 
most complete treatments is found in Clay and Medwin [Ref. 3]. 
The resonance frequencies of air bubbles of given radii at 10 
cm depth of fresh water as calculated using Equation 1 [Ref. 
3:eq. 6.3.11:p. 197] are given in Table 1. Here, b and B® are 
functions that correct for the adiabatic assumption and 
surface tension. These two corrections are counteracting and 
result in very little change to the value calculated by 


ignoring them. Nonetheless, they were included for accuracy. 


1 3YDR pa {3 
pees _- eel (1) 
2a ax 
where, 
a= bubble radius in cm 
Y= ratio of specific heats, C,/C,; taken to be 1.402 
P, = ambient pressure on the bubble corrected for 


depth in dynes/cm 


Pp, = ambient density of the fluid surrounding the 
bubble taken to be 1.03 gm/cm 


see [Ref. 3] for expressions for b and 8B. 


TABLE 1* RESONANCE FREQUENCY FOR A GIVEN BUBBLE RADIUS 
INCLUDING CORRECTIONS FOR DEPTH (10 CM), SURFACE 
TENSION AND THERMAL CONDUCTIVITY 


Bubble Radius, a (pm) Resonance Frequency (kHZ) 


70 
50 
40 
30 
20 i Ze 
10 
g 
8 
7 
6 
5 
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618. 


* The data in this table is included in the Easy Reference 
Table in the Appendix. 


Factors directly effecting the resonance frequency include 
the volume and type of gas, properties of the fluid medium and 
depth of the bubble. For this experiment the gas was 
generated by electrolysis. Various gases can be generated by 
this method depending on the concentration of chemicals, such 
as salts, in the water and the voltage applied. The gas 
produced in this experiment was probably hydrogen, H),. 
Chlorine, Cl, and Oxygen, O, are the other possibilities. 
Since these gases are all diatomic, they share some common 
properties. The principle property of interest is the ratio 
of specific heats, which is listed in Table 2 for these gases 
and air. Also listed is the approximate error induced by 
using the specific heat ratio for air in each case. Since 
these errors are small and many bubble properties have been 
determined and published for air bubbles in water, air was 
assumed to be the gas in the bubbles for all calculations. 
This also allows these results to be more quickly related to 


any fresh water air bubble situation. 


TABLE 2 SPECIFIC HEAT RATIOS FOR SOME GASES 


GAS SPECIFIC HEAT RATIO RESONANCE FREQ. ERROR FROM AIR 
Air 1.402 0.00% 
H, Aa 0.28% 
0, 1.40 0.07% 
1. 1.34 2.24% 


Other factors which effect the resonance frequency, such 
as the depth, pressure, radius and fluid density, are all 
interrelated. Assuming the mass of gas in the bubble stays 
constant, a depth increase raises the pressure and decreases 
the bubbles radius. All of these changes contribute to an 
increase in the resonance frequency. On the other hand a 
fluid density increase lowers the resonance frequency. A gas 
with a larger specific heat ratio would have a higher 
resonance frequency. For this experiment, all resonance 
frequencies are given for air bubbles at 10 cm depth of fresh 
water with a surface pressure of one atmosphere. 

The response of a bubble to an excitation is most dramatic 
when the frequency of the excitation is exactly equal to the 
resonance frequency of the bubble. The bubble's response at 
resonance is much greater than the response off resonance as 
shown in Figure 1 (Ref. 3:Figure 6.4.1: p.204}. 

Due to the large amplitude oscillations of the bubble 
driven at resonance, the higher order terms in the bubble's 
equation of motion are not small with respect to the linear 
terms. These higher order terms are necessary for the dual 


frequency method to work. 


B. DUAL FREQUENCY EXCITATION AND RESPONSE 
References 1 and 2 give extensive treatments of the 


oscillations of a bubble under the influence of two impinging 
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sound fields of different frequencies, f, and f,. f, 1s taken 
to be a high frequency, much higher than the bubble's 
resonance frequency, yet with a corresponding wavelength 
longer than the diameters of the bubbles of interest. f, is 
a lower frequency, at or near the bubble's' resonance 
frequency. As the bubble enters the sound fields, it is 
excited into large amplitude resonant oscillations by f,. It 
is also excited by f,, although at a much lower amplitude. 
The result is that the bubble extracts energy from these sound 
fields and reradiates sound as a monopole source. Because the 
large amplitude oscillations are nonlinear, the bubble 
radiates sound at f, and f,, as well as at the sum (f, + f,) 
and difference (f, - f,) frequencies. All other possible 
combinations are also radiated, but at much lower amplitudes. 
This natural frequency mixing is the cornerstone of the dual 
frequency method. 


1. The Resonant Bubble's Response to Dual Frequency 
Excitation 


Isolating our attention to the sum frequency, the 


pressure radiated by a bubble is given by Equation 2 [Ref. 1]. 


Pi P2 
P+ = ——___—_ (2) 
pawWr* Or 


where, 
p = density of fluid 


Pp, = acoustic pressure amplitude at frequency f, at 
the location of the bubble 


P, = acoustic pressure amplitude at frequency f, at 
the location of the bubble 


p, = acoustic pressure amplitude at the sum frequency 
Cee ate ene hocadelen Of Ene receiver 


W, = resonance (angular) frequency of the bubble 
§ = damping coefficient 
r= distance from bubble to receiver 


It has been assumed that f, is equal to f,, the bubble's 
resonance frequency. The sum frequency pressure amplitude is 
inversely proportional to the distance, r, from the bubble. 
The damping coefficient, 6, is claimed to be almost 
independent of frequency in Reference 1. Values of 6 for 
given bubble sizes at sea level are taken from Reference 3 and 
listed in Table 3. Using Tables 1 and 3, the density of the 
fluid and an assumed distance from the bubble to the receiver, 
the sum frequency pressure amplitude at the receiver can be 
calculated for the bubble sizes of interest in terms of p, and 
p,- This information is shown in Table 4. 

MmepcEEomcuay task similar to that o£ the carrier 


frequency of a radio transmission. Hence, in this thesis f, 


TABLE 3* DAMPING COEFFICIENT FOR SOME BUBBLE SIZES (NUMBERS 
READ FROM FIGURE 6.3.1, P. 199) REVERE Cee 


Bubble Radius, a (pm) § (Damping Coefficient) 
100 0:.0:75 
70 0.085 
50 0.095 
40 0.100 
30 Ozrls 
20 On 13 
10 O16 
<10 Oveal6 


* The data in this table is included in the Easy Reference 
Table in the Appendix. 


10 


TABLE 4* RATIO OF SUM FREQUENCY PRESSURE AT THE RECEIVER FACE 
TO THE PRODUCT OF SOUND FIELD PRESSURES p, AND p, AT 


THE BUBBLE 
pt 
Bubble Radius, a (pm) s(x. :10°’Pa™') 
P,P2 
100 43.2 
Eo ATE 
ae dgjincA 
oe dered 
he 8.79 
a 5.25 
~ Zeaic 
S 1.94 
; eae 
g 1.50 
° meee 
2 1.05 


* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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is called the carrier frequency. It can be any convenient 
frequency high enough with respect to the resonance 
frequencies of the bubbles of interest that the harmonics of 
these bubbles are significantly less than the expected 
difference frequencies. This prevents confusing a harmonic 
with the difference frequency. The bubbles of interest must 
have diameters small with respect to the wavelength associated 
with f,, so that it oscillates uniformly. 

The bubbles of interest to this investigation have 
radii from 100 down to 7 um which corresponds to resonance 
frequencies from approximately 20 to 400 kHz. A reasonable 
minimum wavelength for these bubbles is 400 um, which 
corresponds to a maximum carrier frequency of 3.7 MHz. This 
ensures the acoustic wavelength is at least four times the 
largest bubble radius. Assuming harmonics higher than the 
third or fourth will be insignificant leads to a minimum 
carrier frequency of about 5 times the highest bubble 
resonance or 2.0 MHz. The actual value of f, used in) tie 
experiment was 2.65 MHz which is near the center of this 
frequency range. 

The purpose of f, is to excite the bubbles into Pama 
amplitude resonant oscillations. For this reason, 9 gene 
called the excitation frequency. It must encompass the range 
of resonance frequencies of the bubbles of interest. fThis 


means that f, must cover the full range from 30 to 400 kHz 


ie 


continuously. Two schemes for producing f, are possible. 
One as sweeping the frequency repeatedly over the desired 
range. This method will periodically excite any bubble in the 
range, but only for a short time. Another option is 
broadcasting band limited white noise covering the proper 
frequency range. This more continuously excites any bubble 
in the range at its resonance frequency. This second method 
was chosen so that all bubbles will always find their 
resonance frequency in the f, sound field. This method 
requires more power since it broadcasts over a large bandwidth 
continuously at moderate amplitude. 

Using noise to ensure the resonance frequencies of 
interest are present also ensures that additional frequencies 
are present. This does not confound the dual frequency method 
though, since each bubble naturally responds strongest to its 
resonance frequency. Due to the large quality factors, the 
bubbles will vibrate strongly at their resonance frequencies 
and greatly ignore the sound energy at other nearby 
frequencies. This will allow the detected sum frequency to 
Panpoint the resonance frequency, f,, of the bubble that 
Peeemated 1t. f= £.., - Lcarcier- 

As bubbles get smaller their acoustical cross sections 
(absorption, scattering and extinction) become smaller. To 
keep smaller bubbles excited strongly enough to produce 
observable dual frequency mixing, coloring the noise sound 


field is necessary. A ramp function that causes more sound 


in 6) 


energy at the higher frequencies alleviates this problem 
somewhat. 

A relatively large bubble responds strongly at its 
resonance frequency and also scatters higher frequencies. The 
scattered amplitude can be higher than the amplitude due to 
a bubble resonant at the higher frequency due to the larger 
bubble's large scattering cross section. This leads to an 
ambiguity for other methods as shown by curve b of Figure 2 
(Ref. 1]. For a uniform distribution of bubble radii in the 
mixture, a high bubble resonance frequency may be detected at 
high amplitude due to either, (1) a resonant bubble or (2) a 
much larger scattering bubble. Bubble scattering detection 
techniques looking for the bubble's resonance frequency 
directly cannot be sure which of these mechanisms 1s causing 
the signal without more information. Reference 3 discusses 
this difficulty in more detail. 

The dual frequency method is not susceptible to this 
problem. Mixed frequencies, (f, + f,) and (£, = £5))3ame 
produced only with the large amplitude whole body oscillations 
that occur around resonance. Therefore, since mixed 
frequencies are not simply a scattering phenomenon, it becomes 
unimportant that multiple scattering may be occurring. The 
mixed frequencies are dependent on the bubble's resonance 
frequency, f,, and become (f, + £,) and (f, — £,). Sosa sical 
received at these frequencies indicates a bubble with 


resonance frequency, f,, 1s in the sound field. Even thouaa 
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the scattered signal amplitude peaks, dips and then rises in 
a ramp fashion above the peak as bubble radius increases, as 
shown in curve b of Figure 2, the sum frequency signal is 
strong for only one bubble size, as shown in curve a. 


2. Dual Frequency Interaction for Solids and Non-Resonant 
Bubbles 


Beyer [Ref. 4] has gathered a lot of information 
concerning the nonlinear interaction of dual frequency sound 
fields. In this compilation of results from his own work and 
others, he presents information on the generation of sum and 
difference frequencies by the interaction of the fluid medium 
alone, a solid in the fluid and a bubble in the medium. These 
treatments allow us to estimate a maximum expected signal for 
cases without resonant bubbles. 

a. Sound Interaction with Sound 

In an experiment by Jones and Beyer it was shown 
that no sum or difference frequency 1s generated when two 
sound beams cross at right angles. In fact, sum “amg 
difference frequency signals are generated only when the beams 
are collinear. This case is termed the "parametric end fired 
array." It results ina highly directional sound beam. The 
pressure of this sound beam is proportional to the pressures 
of the primary frequency signals and the frequency separation 
between these signals. To produce the signal, the sound beams 
must be collinear for some distance. As a rule this process 


is less than one percent efficient. 
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b. Sound Interaction with Solids 

The first effect to consider in the presence of 
a solid is scattering. The direction of scatter depends on 
the size, shape and orientation of the solid. To consider the 
parametric end fired array, assume that significant portions 
of both carrier and exciter are scattering toward the 
receiver. Although the receiver is very close to the sample 
volume, some interaction may occur during the interval. The 
resultant signal would have, at most, less than one percent 
of the scattered power. This power estimate is generous 
considering the short interaction length. 

The second effect of solids in a dual sound field 
is their motion due to the radiation forces of the two sound 
fields. The motion of the solid induces the sum and 
difference frequency in the medium it contacts. To have this 
effect, the sound fields must be significantly intense to move 
solids at both the carrier and exciter frequencies. 

c¢. Sound Interaction with Non-Resonant Bubbles 

Beyer gives theory and experimental results for 
an air bubble interacting with dual frequency sound. The 
frequencies used were so high (5 and 7 MHz) that the bubbles 
were not resonant at these frequencies. The theory considers 
the interaction of sound scattered by the bubble. Scattering 
is the primary interaction for bubbles larger than resonance 
Size. For the case presented, a 35 um bubble is ensonified 


by 7 MHz at 2.45 x 10° Pa and 5 MHZ at 3.36 x 10° Pa. These 


AG! 


frequencies are scattered omnidirectionally. The resulting 
pressures 48 cm away are 7 Pa and 6 Pa, respectively. The sum 
Signal arrives with a pressure of 0.07 Pa. 

This is interesting compared to the same analysis 
for a solid sphere. Here the sum frequency pressure is 
proportional to the square of the solid particle's radius and 
varies with direction. The maximum sum frequency signal 
exists directly down range from the sphere on a line with 
either source transducer. For a sphere over 90 times larger 
than the bubble mentioned above and identically ensonified, 
the maximum sum frequency pressure was 0.007 P, at these two 
points. The minimum signal was between these points at a 
pressure of 0.00007 Pa. The placement of the receiver between 
the two sources ensures that we avoid the maximum values of 
the scattered signal pressure. 

d. Overall Analysis of Solid Interference 

Clearly there is a limit to the amount of solid 
material that can pass through the dual frequency sample 
volume without great effect. Much, however, can be said for 
this method's insensitivity to solids. First, the above 
experiments were performed at greater drive pressures than 
required for the dual frequency method. Lowering the drive 
pressures lowers these nonlinear effects. Second, note that 
the solid has much less sum signal than the non-resonant 
bubble. The non-resonant bubble has much less sum signal than 


the resonant bubble, as shown in Figure 2. It follows that 
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a solid particle will have a sum signal that is much less than 
that for a resonant bubble. This provides the means for 
discrimination between bubbles and solids. 

The final argument to separate the resonant 
bubble's sum frequency signal from those of scattering 
phenomena is the frequency selectivity. The excitation signal 
is relatively low power over a broad band. Whereas the solid 
would scatter all frequencies uniformly, the bubble will be 
fundamentally excited by the presence of its resonance 
frequency. The whole body oscillations then mix the BEDmIenis 
resonance frequency with the high power carrier. This sum and 
difference frequency Signal will be easily distinguished from 


the other noise. 


C. ADVANTAGES AND DISADVANTAGES OF THE DUAL FREQUENCY METHOD 
The dual frequency method for bubble detection and 
resonance frequency determination has its good and bad points. 
The following is a quick summary of some of the key issues. 
* Advantages 
- conclusive resonance frequency determination 
- discrimination between solids and bubbles 
* Disadvantages 
- complexity of implementation 


- lower signal to noise levels than direct resonance 
frequency methods 
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* Other Effects 
- rectified diffusion 
- radiation pressure 
1. Advantages 

The foremost advantage of dual frequency bubble size 
determination is the single peaked response shown in Figure 
2. Since this is not a simple scattered sound measurement 
technique, scatterers do not dominate the results. Large 
bubbles respond at their resonant frequency. They do not 
generate the same signals as small bubbles. This advantage 
is the same reason why solids do not give false detections. 

Since this iS an acoustic method uSing actual 
resonance phenomena, true resonance frequencies of skin 
covered bubbles can be determined by experiment. This could 
be important for certain types of plankton or other resonant 
scatterers. The contribution of the skin need not be 
calculated or assumed. 

2. Disadvantages 

One disadvantage is the complexity of the experimental 
set up and support equipment. Function generators must be 
capable of producing colored noise and narrow band high 
frequency. Frequency spectral analyzers must be capable of 
quickly and repeatedly analyzing the received signal over a 
wide range of frequencies. The transducers must be placed 


correctly to overlap the sound fields and receive a signal. 
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The second major disadvantage is the relatively low 
signal to noise level. The amplitudes of the sum and 
difference frequency signals are very small compared to the 
amplitudes of the signals received at the carrier and 
excitation frequencies. The sum and difference frequency 
generation mechanisms are nonlinear interactions and weaker 
than the linear interactions. 

3. Other Effects 

Another problem might be that the sound field could 
change the nature of the bubble being sized. Rectified 
diffusion, a process by which an oscillating bubble actually 
grows, will change the resonance frequency. Then the detected 
resonance frequency would not be the original undisturbed 
frequency. Crum (Ref. 5] gives a detailed discussion of the 
subject of rectified diffusion. The pressure amplitude 
necessary to initiate rectified diffusion is [Ref. 5:eq. 19:p. 


Da? | 


(pa? Wr?) ?[ (1-(W?/WR?) )74+67W? /We® J [1+ (20/ap-)-ca/co] 
Derr) = ee 
(344K) (ca/co)~{[3(n-1) (3n-4) /4 }+(4-3n)K}[1+(20/ap..) J, 


where, 


w= frequency of excitation 


Q 
li 


surface tension of the liquid 


DP. ambient pressure 
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Cc. = concentration of dissolved gas in the fluid far 
from the bubble 


Cc, = equilibrium concentration of gas in the fluid 
K = thermal conductivity of the gas in the bubble 
n = polytropic exponent 
Py, = ##%threshold acoustic pressure at excitation 
frequency. 


Table 5 shows the results of some lengthy computations 
using Equation 3. Comparing the actual pressures used to 
these threshold values verified that rectified diffusion was 
not a problem with this experiment. The other saving point 
is that rectified diffusion takes hundreds of seconds to cause 
an appreciable change in bubble size as shown in [Ref. 5:Fig. 
7:p. 222]. This is much longer than bubbles are expected to 
remain in the sound field. 

Radiation force and streaming are two nonlinear 
effects discussed in Reference 4. These mechanisms have 
effects on the bubbles transiting the dual frequency sound 
field. The radiation force causes the bubbles to move away 
from the regions of high acoustic pressure. Streaming causes 
the bubbles to move away from the source transducers when in 
the sample volume. These motions make it difficult to 
maintain a small bubble in a small sample volume long enough 
to detect it. The effect of both of these mechanisms can be 
reduced by lowering the carrier and excitation pressure 


levels. Having a large sample volume with respect to the 


ee 


TABLE 5. THRESHOLD PRESSURES FOR RECTIFIED DIFFUSION OF 
VARIOUS AIR BUBBLES IN AIR SATURATED WATER 


(eye, = 2) 
Bubble Radius Resonance Threshold Frequency 
(sm) Frequency (KHZ) Pressure (Pa) (KHZ) 
100 31.4 1438 31.4% 
50 Ge. 255) 62, 1* 
10 SUS 8676 303.8% 
5 618.0 1515000 2650.0** 


* W/W, = 1 


*k = Carrier frequency 
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distance a bubble moves under the influence of these forces 


would also mitigate their effect. 


D. BUBBLE RISE TIME SIZING 

In order to verify that a bubble's acoustically indicated 
size is correct, some other method of determining the bubble 
Size is necessary. A relatively easy, straightforward and 
accurate method is to measure the rise time of the bubbles as 
they rise through a known distance in still water under the 
influence of buoyancy. Reference 6 has documented this 
technique. 

Essentially, a bubble rises at constant speed due to the 
balance between its buoyancy and its drag. These forces have 
different dependencies on radius. Hence, the rise speed of 


bubbles varies with radius; the larger the bubble, the faster 





it rises. The terminal rise speed of a bubble is given by 
(Ref. 6] 
8 ag }3 
pies (4) 
SuCy 
where, 
U= rise speed 
g= acceleration due to gravity 
C, = drag coefficient 
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The values determined for time to rise one inch for 
various bubble sizes are shown in Table 6. The Schiller and 
Nauman drag law was used for these computations as given by 


Equation 5 [Ref. 6:eq. 4:p. 6]. 


24 
Gowen Reo oe 7) (2) 
Re 


where, 
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Reynold's number = 2aU/V, 


V,; = kinematic shear viscosity of fluid 


In practice, there are a few difficulties with this 
method. The experimenter must be able to monitor the moving 
bubble. Thirty micron radius bubbles are about the smallest 
easily timed, with good lighting and clear water. There is 
some parallax since the scale cannot be extremely close to the 
rising bubble. With large bubbles, which move fast, observer 
response time is a factor. The experimenter cannot be 100 
percent certain that the bubble timed is not accompanied by 
other bubbles that effect the acoustic output. 

Overall, the rise time method is helpful. rt 
allows at least a rough comparison of the acoustically 


determined resonance frequency with an expected value for 


ZO 


TABLE 6* TIME TO RISE ONE INCH FOR VARIOUS BUBBLE SiZee 


Bubble Radius, a (pm) Time to Rise One Inch (sec) 

100 1.44 
70 PADD), 
50 4.61 
40 64°95 
30 i240 
20 26.5 
18, AsO 

9 129 

8 163 

7 213 

5 416 


* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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those bubbles that can be timed. Even when it cannot be used 
for every bubble, it gives the experimenter confidence in the 


dual frequency method. 
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III. PROBLEM APPROACH 


The goal of this thesis is to develop and test a system 
which uses the qual frequency method to detect and size 
bubbles. In order to accomplish this goal, the following 
capabilities had to be incorporated into the experimental 
system. 


* A method of generating bubbles and placing them in the 
sample volume. 


*x* A method of placing solid particles in the sample volume. 


* A method of producing a sound field of sufficient pressure 
amplitude. 


* A method of receiving the pressure signals from the bubble 
and determining the frequency spectrum of that signal. 


These methods are discussed in this chapter. 


A. BUBBLE GENERATION 

A method for generating bubbles of very small size is 
necessary for this study. Two basic generation methods were 
tried. In the first, bubbles were generated by pressurizing 
a pipette, made from a glass rod with a very narrow axial 
hole. The pipettes were very delicate, produced an audible 
sound during bubble generation and produced bubbles much 
larger than desired for this study. For these reasons, this 
method was abandoned and an electrolytic method was developed. 


Figure 3 shows the basic bubble production setup. 
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Hewlett Packard 
Tap 3314A 
Water a Function Generator 





Figure 3 Electrolytic Bubble Generation Setup 


Using the arbitrary (ARB) wave function mode on the 
Hewlett Packard (HP) 3314A, a negative electrical pulse was 
generated. The pulse was applied to a small wire, all but the 
very tip of which was insulated. A return wire with a bare 
coil tip completed the circuit back to the function generator. 
Ordinary (and by no means pure) tap water served as the 
electrolyte between the two leads. Very small bubbles would 
form and rise from the small wire during the negative pulses. 

The nature of the bubbles was easily varied, but not 
controllable. Because the small wire tip corroded and 
collected deposits rapidly with use, the bubble size and 
quantity would vary with time, even with all of the parameters 
held constant. The use of a non-corroding wire, platinum for 
instance, might have solved this difficulty. 

The number of possible variations offered by this setup 
proved useful. By adjusting the time scale on the ARB cycle 


the bubbles could be made to appear in quickly recurring 
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bursts having small numbers of small bubbles or in groups 
having a large number of larger bubbles with a larger interval 
between groups. As a rule to obtain smaller bubbles: (1) use 
a thinner wire, (2) use a shorter duration negative pulse, (3) 
use a lower voltage negative pulse. Conversely, larger 
bubbles were more likely with thicker wire, higher voltages 
and longer duration pulses. 

It was impossible to generate only one size bubble over 
any length of time with this method. With several small 
bubbles present, two would often join to form a larger bubble. 
The larger bubble would then rise faster, overtake and join 
another bubble. To prevent this, an on-off switch was placed 
in the circuit. The current was interrupted after a bubble 
was formed to limit interference with other bubbles. 

1. Variations of Bubble Control 

a. Grid 

Since the position of the rising bubbles is 
important to this experiment, some way is desired to move the 
generation point a known amount. One solution is a grid of 
32 wires terminating on a small blank circuit board. Arranged 
in 4 columns and 8 rows with the tip of each wire about one 
millimeter from its neighbor, the source of bubbles could be 
moved by switching one wire off and another on. Also a 
massive cloud of bubbles could be achieved by turning on all 


or some of the wires at the same time. 
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b. Bubble Size Discrimination Plate 

In order to narrow down the size range of bubbles 
which passed through the sound field, a small horizontal water 
jet was produced near the bubble generator. The jet was 
established by connecting a water reservoir to a small hollow 
glass tube. The speed of the jet could be controlled through 
the height of the reservoir. As bubbles were generated they 
were swept downstream by the jet. The length of travel was 
inversely proportional to the size of the bubble. Eventually 
the bubbles would rise out of the jet and continue up to Sie 
surface. A plastic plate with one small hole was placed above 
the jet. This arrangement allowed only those bubbles leaving 
the jet at the position of the hole to rise through the sound 
field, as shown in Figure 4. Although helpful in narrowing 


down the range of bubble sizes, some variation was still 


apparent. 
Plastic Plate 
Ties Sie 
© 
© 2 O 
. 2 © Tube To Water Tank 
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Water Jet O35 
© 
———— TO HP 3314A 
Wire 


Figure 4 Bubble Size Discriminator Plate 
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c. Blow Pipe 

There are a number of difficulties associated with 
getting very small bubbles into the sound field. Worst of 
all, they are just about impossible to see. They drift away 
from the acoustic sample volume very easily, and they rise at 
very slow speeds, several minutes per inch. In order to get 
them into the sample space, they had to be convected. To this 
end, the bubbling wire was inserted into a narrow glass tube. 
The tube was pressurized by a water reservoir with about one 
foot of static head. The water was allowed to flow 
continuously while the bubbles came in short pulses. The tube 
could be (relatively easily) positioned so that the convected 


bubbles would pass through the sample volume. 


B. SOLID INJECTION 

In order to show that the frequency mixing due to solid 
particulates in the system is dramatically different from that 
due to bubbles, and that therefore the system would show the 
ability to discriminate between a bubble and a solid, such 
solid particles had to be introduced into the sample volume. 

Three basic types of solid particles were injected. The 
first, dry sand, was dropped directly over the dual sound 
field via a dry funnel and allowed to fall to the tank bottom. 
It was never certain whether the particles fell directly 
through the most intense portion of the sound field or not, 


due to their long fall and fluttering motion. 
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Wet sand was injected in a water slurry via a long wet 
funnel terminating just above the dual sound field. Wetting 
the sand prior to injection reduced the number of visible 
bubbles entrained with the sand. 


Thin wall hollow glass beads were also injected into the 


sound field. Due to their extreme buoyancy they had to be 
convected. A glass tube was aimed at the dual sound field 
horizontally and pressurized with a bead slurry. The beads 


would pass into the sound field horizontally and rise to the 
surface. The glass beads ranged in size from about 5 to 50 
um radius, as determined by microphotography. The sand 


particles ranged from 1 um to 1 mm in largest dimension. 


C. TRANSDUCER CHARACTERISTICS AND CONFIGURATION 

Due to the range of frequencies used in this experiment 
several types of transducers were needed. The lower frequency 
driving transducer, called the exciter, needed significantly 
different parameters than the higher frequency driving 
transducer, called the carrier. Since the sum and difference 
frequencies were close to the carrier frequency, the receiving 
transducer, or receiver, was very similar to the carrier 
transducer. Other transducers were used for special purposes, 
such as calibration, but were not part of the dual frequency 
system. This section will present each transducer's 
characteristics and the methods used to determine those 


Siharacteristics. 
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1. Types of Transducers 

Two 6 cm diameter focused piezoelectric transducers 
were used for the carrier and receiver. They were chosen for 
the acoustic gain advantage of focusing. The acoustic 
pressure at the focus iS many times that at the face of the 
transducer. The serial numbers of these transducers are C5575 
and C5574. They were manufactured by Harisonic. Nominal 
specifications include a 2.25 MHz resonance frequency and a 
3 inch focal length. 

A Panametrics V30u 3 cm diameter cirewiars 
piezoelectric disk transducer was used as the exciter. This 
choice was made since it had good frequency response over the 
range of 30 to 400 kHz. This transducer had a source strength 
that mainly increased with frequency over this range. That 
was ideal for this purpose. The serial number of this 
transducer was 93598. To assist in reciprocity calibration 
a Similar transducer, serial number 93596, was used. 

To determine the transducer beam patterns and assist 
in correctly positioning the transducers in the system, a very 
small hydrophone, made by Specialty Engineering, was used. 
The hydrophone looks like the tip of a very sharp pencil, with 
all dimensions of the ceramic receiver less than a millimeter. 
A preamplifier supplies enough signal boost to give 
significant oscilloscope readings. This hydrophone was used 
to probe the sound fields for relative measurements. No 


attempt was made to calibrate the probe. 
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Finally, an EDO corporation model 6600 spherical 
piezoelectric hydrophone was used for the reciprocity 
calibration of the exciter. 

2. Beam Patterns 

The beam patterns of the exciter, carrier and receiver 
are important since their intersection determines the size of 
the sample volume. Also, an accurate knowledge of the beam 
pattern is crucial to the calibration of the focused 
transducers. The methods used to measure the beam patterns 
of the transducers are described in this section. 

a. Beam Pattern Method 

In order to obtain an accurate beam pattern, the 
transducer and probe hydrophone were mounted to a common frame 
via a network of micrometer positioners. Thus, the probe 
could be moved known distances, relative to the transducer, 
in small steps. The setup allowed relative motion in the 
three orthogonal directions of vertical, range and cross 
range. 

The first step in measuring the beam patterns of 
the focused transducers was to locate the position of the 
mOCcuUs . Then readings of distance and probe output voltage 
were taken relative to the focus. At any given frequency, the 
probe output voltage was assumed to be proportional to the 
acoustic pressure at the probe. Normalized pressure values 
were obtained by dividing each value in a data set by the 


maximum value in that data set. 


eile) 


The exciter beam pattern was much broader than 
those of the focused transducers. Cross range data for the 
exciter was taken at a few ranges of interest. Since the 
exciter's beam pattern is, theoretically, narrower at higher 
frequencies, the data was taken at 400 kHz. This verified 
proper coverage at the highest exciter frequency of interest. 
The patterns at lower frequencies were then guaranteed to be 
sufficiently broad. 

b. Beam Pattern Results 

At 400 kHz the exciter had a very smooth beam 
pattern. At a distance of 4.0 cm the 3 dB beam width was 
about 1.0 cm. Figure 5 shows the normalized pressure versus 
cross range at a range of 4.0 cm. The variation of pressure 
with range was minimal for ranges 4 cm to 6 cm. At a range 
of 10 cm the maximum pressure was approximately 3 dB less than 
that at 4.0 cm range. Figure 6 shows the area within which 
the pressure waS greater than the 3 @B less than the maximum 
value (0.71 p,,,). AS will be shown below, this beam pattern 
is sufficient to ensonify the intersection volume of the 
carrier and receiver, which limited the size of the sample 
volume. 

The beam patterns of the focused transducers were 
much more complicated. The acoustic pressure changed rapidly 
with distance, especially at the focus. Figure 7 shows the 


variation of normalized pressure with range along the axis for 


36 


Normalized Pressure 


0.8 


0.6 


0.4 


2 


10 20 S10. 


Cross Range (mm) 


Figure > Exciter Cross Range Pressure 
Variation at a Range of 4.0 cm 


a) 


25 


20 


Cross Range (mm) 
UI 


4 6 8 10 


Range From Exciter Face (cm) 


Figure 6 View of Exciter Main Beam. Area inside curvesmie 
High Pressure Region. 


38 


Normalized Pressure 


60 70 80 
Range (mm) (Relative Scale) 


migure / Focused lransducer s Pattern of Pressure versus 
Range at the Focus 


39 





90 


transducer C5574. Both focused transducers were similar in 
this regard. At the focus, the pressure stays at or near the 
maximum for about 5 mm of range, and at least half of the peak 
value for 10 mm. Figures 8 and 9 show the vertical variations 
of normalized acoustic pressure at the focus of C5574 and 
C5575, respectively. The pressure stays at or near its 
maximum value for only 1 mm vertically. Comparing these 
graphs for the two transducers reveals that C5575 focuses a 
little tighter than C5574. The cross range beam patterns are 
very Similar to the vertical patterns and are not shown. 
Although the beam patterns are axially symmetric in general, 
there are a number of small variations in the pressure field. 

Figures 10 and 11 show the horizontal plane (cross 
range versus range) containing the focal axis for C5574 and 
C5575, respectively. Isobaric lines surround areas within 
which the pressure is at least as great as the isobaric 
pressure. This is just another way to show that pressure 
varies ten times faster with cross range than with range. To 
compare the graphs, note that although the isobaric lines are 
labeled differently in units of dB relative the maximum value 
for that plot, the inner oval for C5574 represents the same 
pressure as the inner oval for C5575. These plots were 
instrumental in understanding the focused transducers' sound 
fields. They provided key information for the calibration 


processes as shown in the next section. Further detailed 
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analysis of the resulting sample volumes is given after the 
calibration section. 
3. Transducer Calibration 
The beam patterns in the previous section were all 
plotted relative to some maximum value for that transducer. 
To apply the dual frequency theory and calculate the expected 
sum and difference frequency signal pressures, the carrier and 
exciter pressures must be determined. Also the receiver 
sensitivity must be known. To this end, several calibration 
methods were used. This section presents a brief outline of 
each method and the results obtained. 
a. Methods of Calibration 
(1) Reciprocity. References 7 and 8 present a 
significant amount of information on reciprocity calibration 
of transducers. If done correctly, this type of calibration 
is a prime standard of transducer performance. The key to the 
reciprocity method is two fold, as Reference 9 points out. 
First, for a linear, passive and reciprocal transducer there 
is a simple relationship between the transducer's receiving 
sensitivity, M, and its transmitting response or source 


strength, S; 


yee (6) 
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where, 
M has units of voltage/pressure 
S has units of pressure/current 


J has units of (voltage) ‘(current) /(pressure)*. 


J is called the reciprocity parameter and depends on the 
nature of the calibration method. Several methods of 
reciprocity calibration are presented in Reference 8 with 
their associated reciprocity parameter. The second key to the 
reciprocity method is that the beam pattern of a transducer 
is the same whether it is acting as a source or receiver. 
This property was assumed in order to determine the beam 
pattern of the receiver in the last section. 

The full reciprocity method requires one 
transmitter, T, one hydrophone, H, and one reversible 
transducer, R. Three basic steps exist in the experiment. 
First, place the transmitter, T, a known distance, d, from H. 
Record the open circuit voltage response from H, V,,, for a 
known current, i, into T Second, replace H with R. Record the 
open circuit voltage out of R, V,p;, for the same current into 
dus Third, replace T with R and return H to its original 
position. Record the open circuit voltage out of H, Vip, for 
a known current, ip, into R The hydrophone sensitivity, ' M,, 


can be calculated according to 
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ae = ee (7) 
Vrr lr 
where, 
J = 2id/p,c, for spherical spreading 
A= acoustic wavelength 
p,.c = characteristic impedance of the medium. 


A transducer's source strength, S, is a measure of how much 
pressure, p, it produces at a distance of one meter for a 


given current, i. When p is measured at a distance, d, 





(8) 


Sensitivity, M, is a measure of the open circuit voltage, V, 


produced for a given pressure, p, at the transducer face, 


< 


gol | 


For this reciprocity calibration it is assumed that the 
hydrophone is small enough to not disturb the sound pressure 
field. Thus, the pressure in the above two equations can be 


equated, giving 
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V d 
S = 
Mi BRC) 


(10) 





This equation is used to calculate the exciter's source 
strength after the EDO hydrophone's sensitivity is determined. 

(2) Two Transducer Reciprocity. A short cut can 
be made if two identical reciprocal transducers are to be 
calibrated. Here, only one transducer arrangement is 
required. The transducers are placed a known distance, d, 
apart. Then the open circuit voltage of one is measured while 
a known current is put through the other. The sensitivity and 
source strength of both transducers are given by Equations 11 


and 12. 


<j 
Nim 
= 


M= = Sanco a= (11 and 12) 


# 
Cy | 


The disadvantage of this method is that the transducers must 
be exactly identical in sensitivity and source strength. This 
must be shown by another method before the two transducer 
reciprocity becomes valid. 

(3) Self Reciprocity Method. A third variation 
on the reciprocity method involves only one transducer. This 
method is very similar to two transducer reciprocity, except 
that there is no longer the need for identical transducers. 


The transducer is aimed at a perfect reflector, such as an 
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air-water interface. A pulse of known current is put through 
the transducer. The acoustic pulse is reflected from the 
interface and returns to the transducer. The open circuit 
voltage amplitude of the received pulse is recorded. Care 
must be taken to ensure the transducer is aimed correctly to 
receive the reflected pulse. The source strength and 
sensitivity equations are identical to those for the two 
transducer methods. Since the transducer remains connected 
to the pulse generator while measuring the return voltage, a 
Giode or other high impedance arrangement must be caeaweee 
isolate the influence of the pulse generator's low impedance 
on the open circuit voltage measurement. 

(4) Radiation Force Target Deflection. In this 
method, a perfectly reflecting target is hung in the 
transducer's sound field. The target is deflected due to the 
radiation force. The amount of deflection is determined by 
the’ balance of the radiation force and the weight of the 
target. As shown in Reference 10, the total acoustic power, 
I], acting on the target is determined by the deflection 
distance, e, multiplied by a target specific factor, as shown 


in Equation 13. 


I] = — e (13) 
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where, 


by Equation 


mass of the target 

gravitational acceleration 

speed of sound 

length of target suspension 

deflection distance. 

At the focus, the total acoustic power is given 


14. 


II = | I- dAr (14) 


intensity of the sound over the infinitesimal 
area, dA; 


infinitesimal area in a plane intersecting the 
focus and perpendicular to the transducer's axis 


The intensity is related to pressure, for plane or spherical 


waves, by Equation 15. 


p? 
PoC 
Py 
I] = _—_— dA, (TGs) 
oc 
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Now, define s(x) as the pressure's radial 
shape function which equals the normalized pressure at a 
vertical or cross range distance, x, from the focus. In other 
words, this shape function looks exactly like Figure 8 for 
C5574 or Figure 9 for C5575. The pressure at a radial 
distance, x, in this focal plane, p,(x), can be written in 
terms of the maximum pressure and the radial shape function 


as in Equation 17. 


Pr(X) = PemaxS(X) (er?) 
substitution gives, 

Prmax 09 
I] = 27 S*(x) xdx (18) 

PaC O 

Prmax 
I] = ————— Ave (19) 

pC 

A.;, 1s the effective beam area at the focus. It is obtained 


by numerically integrating the shape function squared over the 
focal plane. 
Substituting the target deflection parameters 


for fl and solving for maximum focal pressure yields Equation 


20. 


Nip 


PoC* mge 
a= (20) 





Prmax = 


ll 


In order to compare this result to that of 
the self reciprocity method, the pressure at the transducer 
face (determined from the source strength) must be converted 
to the pressure at the focus. The conversion is made assuming 
the acoustic energy at the face equals that at the focus. 
This assumption yields 


Dr“Ar 5 
(21) 





Drmax = 
Ae ff 


where, 


A, area of the transducer face 


P, pressure at the transducer face. 


b. Calibration Results 

Some of the methods described in the previous 
section were more suitable than others depending on the 
particular transducer to be calibrated. For instance, the 
full reciprocity method was ill suited for the focused 
transducers because of the rapidly varying pressure pattern 
in the focused sound field and the high frequency. Under 
these circumstance, the assumption that the hydrophone does 
not disturb the pressure field is not valid. The other tree 
methods were viable alternatives. For the exciter, the Self 
Reciprocity and Radiation Force Target Deflection methods were 


not used because of the low pressure amplitudes. 


SZ 


(1) Exciter Calibration. A full reciprocity 
calibration was performed on the exciter, using a second 
transducer of the same design (twin) and an EDO hydrophone. 
This type of calibration results in a value for the 
hydrophone's receiving sensitivity. The exciter's source 
strength can then be calculated by Equation 22. 

Veno d 


Ss = ————— (22) 
Mevo le 1(m) 





where, 
S. = exciter's source strength 
Veno open circuit voltage of EDO hydrophone 
Moy7 =  %§EDO hydrophone's receiving sensitivity 
i, = current input to exciter 
a= distance between source and hydrophone 
The calibration was performed twice at each 
frequency to get two values for the EDO's sensitivity. To 


make the two values somewhat independent the roles of 
projector and reciprocal transducer were exchanged between the 
exciter and its twin. The values for the EDO's sensitivity 
were averaged and compared with expected values provided from 
the manufacturer. The average EDO sensitivity was used in all 
computations. 

The source strengths of the exciter and its 


twin were calculated using Equation 22. At some frequencies 
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the source strengths of the exciter and its twin were close 
to identical. The two transducer reciprocity calibration was 
performed twice in that case. Again, the roles of source and 
receiver were switched to provide some small degree of 
independence between the two runs. For this method the source 


strength is given by Equation 23. 


S = — (239 


where, 
V = open circuit voltage of the receiver 
i = input current to the source 
J = reciprocity parameter, 2,Ad/p,c 


X = acoustic wavelength 


For these calibrations, the received open 
circuit voltage was measured directly on an oscilloscope 
having a 1 megaohm input impedance. The input current was 
determined by measuring the voltage across a 10.1 ohm 
resistor, placed in series with the driving transducer, on 
another channel of the same oscilloscope. The drive signal 
was obtained from a Hewlett Packard (HP) 3314A Function 
Generator. 

The full reciprocity method seemed to be the 
most reliable of these calibrations. The full reciprocity 


calibration results for the exciter were used in all further 


54 


equations requiring exciter source strength. Since the 
exciter was used for the full continuum of frequencies from 
30 to 400 kHz, it was calibrated at 50 kHz intervals from 50 
to 400 kHz. Table 7 contains the results of these 
calibrations. 

The true matter of interest is the pressure 
in the sample volume, which is obtained from Equation 8 
rearranged as p =s i‘1l(m)/r. The drive current is determined 
by measuring the voltage across the 10.1 ohm resistor .as 
discussed earlier. However, the drive signal is broadband 
noise. The oscilloscope cannot be used to determine voltage 
at a chosen frequency in this case. The voltage (rms) was 
measured with the HP 3585A Spectrum Analyzer. The peak 


current is therefore given by Equation 24. 


1 = ——___ (24) 


Using this expression the pressure is found to be, 


Ve Veme S°1i(m) 
ae (25) 
R r ’ 





SYS, 


TABLE 7 EXCITER CALIBRATION RESULTS. 
SOURCE STRENGTH IN UNITS OF Pa/mA 


PULL RECEPROGT iw TWO TRANSDUCER METHOD 


FREO. (KHZ) StpACT ER) S (TWIN) 5 (EXCITE RG S (TWIN) 
50 310 22.6 a <= 
100 1440 1160 eae <a 
150 1860 25,0 =o a 
200 3340 3200 3380 3440 
20 5160 5300 5400 5460 
300 11100 9260 = =a 
Bore 11300 10300 11200 11200 
400 13600 13700 14200 14200 
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where, 


Ving = VOltage determined by spectrum analyzer 


resistance of series resistor 


v8) 
I 


r = distance from the transducer face to the sample 
volume 

Table 8 gives the values of peak pressure for 
the maximum noise signal used to drive the exciter during the 
experiments. It turns out that the exciter has some 
significant variation in output pressure over this frequency 
range. Interpolated values were used for frequencies of 

interest that were not part of the calibration. 

(2) Focused Transducer Calibration. Thestirst 
method attempted in the focused transducer calibration was the 
two transducer method. The position of both transmit and 
receive focused transducers were adjusted to maximize the 
return signal. This position turned out to be such that the 
transducers faced each other with the acoustic’ foci 
overlapping at half the distance between them, as shown in 
Figure 12. The two transducer method was performed twice. 
Switching the roles of source and receiver between C5575 and 
C5574. The obvious problem with this method was that it gave 
no indication of differences between the transducers. The 
result is a single number for source strength and another for 


sensitivity. It was suspected and later proven that the 


a7 


TABLE 8* PEAK EXCITER PRESSURE AT 4 CM AT VARIOUS 
FREQUENCIES WHEN DRIVEN WITH MAXIMUM NOISE 


SIGNAL 
Frequenc KHz Peak Pressure at 4 cm (Pa) 
S)0) LYS 
100 Zo. 
150 244 
200 wae 
250 Bales 
300 980 
350 1000 
400 de. 9 0 


* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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Figure 12 Two Transducer Reciprocity Calibration Using 
Focused Transducers 








2, 


transducers were not identical. This calibration resulted in 
a ball park verification of other methods used. 

Since the impedance of the focused 
transducers was not large compared to 10.1 ohms, and external 
resistor was not used in the determination of input current. 
Instead the magnitude of the free field impedance, |z|, of the 
source transducer was measured separately at the frequencies 
of interest using a HP 4192 LF Impedance Analyzer. With this 
value known, the input current could be deduced from the drive 
voltage. The receiving transducer's open circuit voltage was 
measured directly on the oscilloscope. 

As stated previously, one key to any 
reciprocity calibration is the reciprocity factor, J. For 
small, omnidirectional sources in open water this factor is 


normally given by Equation 26. 


J. = (26) 


where, 
A = acoustic wavelength 


p.c = characteristic impedance of the medium 


The distance factor, dad, 1s included to 
account for losses due to spherical spreading. However, for 
the focused transducers, spherical spreading is not the case. 


The sound field converges to and diverges from a focal point. 
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On the average, this can be equated to a plane wave phenomena 
in the sense that the receiver intercepts essentially all of 
the acoustic energy produced by the transmitter. Reference 
8 gives a brief presentation of plane wave reciprocity. 
Reference 11 provided a verification of the plane wave 
reciprocity factor given in Reference 8 and Equation 27. 

2rd 


Jp = —— (27) 
Fac 


where, 
A = area of transducer face (curved area for focused 
transducers) 
In order to avoid the need for identical 
transducers, the self reciprocity method was used. The 


transducer was positioned vertically underwater to achieve a 
maximum surface reflected signal. This placed the focal point 
at the air water interface as shown in Figure 13. 

A pulse of about 200 cycles at the desired 
frequency in the range of 2.25 MHz to 3.10 MHz was generated 
by the HP3314A every 10 milliseconds. The drive pulse was 
powerful enough to pass through the diode network and drive 
the transducer under test. However, the diodes caused some 
Minor Signal distortion. To account for this, the drive 
Signal was monitored on a spectrum analyzer at each drive 
frequency used. A ratio of signal voltage at the drive 


frequency aS measured on the spectrum analyzer to the waveform 
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Figure 13 Self Reciprocity Method for Focused Transducers 
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peak to peak voltage as measured on the oscilloscope was 
calculated for each drive frequency. This ratio was used as 
a correction factor to the drive voltage measured on the 
oscilloscope. 

The diodes essentially eliminated the HP 
3314A's low impedance from the receiving circuitry. Because 
of losses to the diode circuit, the received voltage measured 
on the oscilloscope with the diode circuit installed was 
noticed to be about 99 percent of the open circuit voltage. 
A factor of 1.01 was used to correct the oscilloscope readings 
to open circuit voltage readings. Since the received signal 
was spectrally dominated by the desired frequency, a second 
frequency correction factor was not required. 

Table 9 shows the results oof this 
calibration. In general, C5575 was a more powerful source and 
a more sensitive receiver than C5574. Notice that the two 
transducer method result of 0.135 mV/Pa sensitivity at 2.65 
MHz is not far from an average of the self reciprocity method 
results for the same frequency. 

Note also the unexpected difference in the 
characteristic rolloff from the nominal resonance frequency 
OG 2.25 MHz. C5575 shows an expected pattern of falling 
sensitivity with increasing frequency beyond this point. 
While C5574 stays virtually flat over the sum and difference 
frequency range. For this reason, C5574 was chosen to be the 


receiver. C5575 provided the carrier signal. 
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TABLE 9 SELF RECIPROCITY CALIBRATION RESULTS FOR VT 
FOCUSED TRANSDUCERS 


C5575 C5574 

* * SOURCE ** SOURCE 
Jeg eo). SENS PTI iy STRENGTH *SENS aL Ty STRENGTH 
(MHz ) (mV/Pa) (Pa/mA) (mV/Pa) (Pa/mA) 
Zeao 0.244 56.3 O209s2 232 
Page 8) OO. 238 5i/ ac0 O20969 23 
Pes SD U.256 56.6 On 09S) 226 
2.40 Uae DO > O203 ou 22am) 
2.45 One si2 D5.) Oh ONS) ae. Loans 
ZO Oizo 0 DosU 0.0936 22.4 
Py esS 9, 0.230 Say! 0.0932 Pa Pe 
200) Oe229 54.8 0.0934 2225 
ZO Or227 54.4 0. 0932 2228 
227 0 0.226 54.1 0709377 22.4 
PRE aS. 0.224 2350 0.0948 OXON 
2.80 Os221 22729 0.0947 2210 
Fi oe. On-2 17 D260 0.0967 2 eel 
Page) 9, O.213 SOAS, O03 66 2a wel 
PA So: 0.207 49.4 Or Oe 2. Suse 
3.00 0.196 46.9 0.0974 230s 
a0 5 0.186 44.5 02-0972 2362 
3210 0.176 42.1 0.0966 23 


* This data is included in the Easy Reference Table in the 
Appendix. 


*k at the transducer face 
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Of particular interest is the pressure of the 
carrier signal in the sample space. The maximum pressure at 
the focus, Pray, is calculated from Equation 28 after 


determining the input current, i, and the effective acoustic 


area, Ag. 


Nie 


(Sa .6serz X i)7 Ar 








Pemax = (28) 
Aer? 
y 
wa PP (5) 
2|Z|l2.6smerz 
where, 
Vp = peak to peak voltage of drive wave form on the 
oscilloscope 
\Zlo.¢smyz = magnitude of free field impedance of C5575 at 2.65 
MHz 
So 65muz = source strength of C5575 at 2.65 MHz 


Since the carrier is monofrequency, the only 
Peeameter that varies is V,,. Table 10 shows some calculations 
for Prewax based on some typical input voltages. The following 
values were used for this table. Sogsmy, = 54. Pa/mA, A, = 


3094.5 mm, Ay = 0.6811 mm, |z]>¢syyz = 135.3 ohms. 
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TABLE 10* MAXIMUM PEAK FOCAL PRESSURE OF THE CARRIER AT 2.65 
MHZ FOR VARIOUS INPUT VOLTAGES AS DETERMINED BY 
THE SELF RECIPROCITY CALIBRATION 


Voltage (Peak to Peak) Maximum Peak Focal Pressure (Pa) 
O..052Z 707 
0.165 2242 
0.41 Bo) 
OES 5 7472 
Bee 8 14940 
2206 28260 
PAPE SS. 34640 
4.05 55020 
5.2 70650 
Sed 111400 

DOs SD 140600 
1a 175300 


* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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Figure 14 Radiation Force Target Deflection Calibration Setup 
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An additional calibration using the radiation 
force target deflection method was performed in order to 
verify the reciprocity calibration results. This test was 
easily performed using the setup shown in Figure 14. The 
target was positioned at the approximate location of the focus 
in order to intercept the entire signal. The deflection was 
measured with a telescope mounted to a micrometer positioner. 
The HP 8640B Signal Generator created the 2.65 MHz signal with 
very fine voltage control and spectral purity. Using Equation 
20, Table 11 was generated using the deflection distances and 
input voltage data from this experiment. This determination 
of peak carrier pressure at the focus compares favorably with 
the self reciprocity method. 

4. Test Setup and Sample Volume 

Due to the highly focused beam patterns, precise 
alignment of the transducers was required. For this reason, 
the transducers were mounted to micrometer positioners which 
allowed adjustments of range, cross range and height. The 
alignment procedure consisted of first placing the probe 
hydrophone in the focal region of the carrier. The receiver 
was then carefully positioned so that its focus also coincided 
with the position of the probe. Similarly, the exciter was 
aimed at the probe to provide maximum excitation signal to the 
sample volume. The bubble generator was then positioned such 


that the bubbles impinged on the probe and thus would go 
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TABLE 11* MAXIMUM PEAK FOCAL PRESSURE OF THE CARRIER AT 2.65 
MHZ FOR VARIOUS INPUT VOLTAGES AS DETERMINED BY 
THE RADIATION FORCE TARGET DEFLECTION CALIBRATION 


Voltages (Peak to Peak) Maximum Peak Focal Pressure (Pa) 


4.0 47080 
5. 0 66580 
6.0 81540 
in W) 88080 
or 105300 
9.0 ¥T5500 
10.0 ro 2 00 
Jobs 10) 148900 
r2.0 166400 
oO 188300 
14.0 Ibe 1 00) ©) 
Po. 0 Zola 7 O00 
iG. © 230600 
20.0 264200 
Zo So9500 
S018) 410400 


* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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through the sample volume. The probe was then removed and the 
experiment proceeded. 

Figure 15 is a full size drawing of a focused 
transducer and its focal region. This representation 
emphasizes how small the focal region is for the carrier and 
receiver. Figure 16 is a similar figure for the exciter. 
Figure 17 shows the three beam patterns overlayed (not to 
scale). The intersection of these three regions is the sample 
volume. The receiver and carrier can be aligned no closer 
than 55 degrees due to the size of the transducers. The 
exciter provides full coverage when properly aimed from almost 
any angle. A value for the sample volume size can be 
determined for a given minimum acoustic pressure assuming the 
receiver and carrier foci overlap. The region can be 
approximately pictured in three dimensions as the intersection 
of two circular cylinders at an angle of 55 degrees. The 
radius of the cylinders depends on the acoustic pressure 
required for the dual frequency method to be effective. In 
the horizontal plane containing the focus, the area, Ay,,, 
through which a bubble must rise to be in the sample volume 


is given by 
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Figure 17 Beam Pattern Intersection 
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where, 


r. = effective radius of carrier cylinder 
re = effective radius of receiver cylinder 
@ = beam pattern intersection angle 


The maximum vertical dimension of the sample volume is the 
minimum of the diameters of the intersecting cylinders, called 
Huax ° 

These equations were used to calculate the values 
shown in Table 12. The sample volume is not isobaric. The 
pressure is maximum in the center. Volumes of decreasing 
pressure surround the focus much like the skin of an onion. 
The geometry gets much more complicated farther axially from 
the focus, especially if the foci are not at the same point. 

Using the known sizes and rise velocities of the 
bubbles, it is possible to determine the length of time a 
bubble is in the carrier field at a given pressure. It is 
also important to realize that it is very unlikely that the 
bubble will pass directly through the focus. Figure 18 is a 
Simplified view of the sample volume cut into volumes of 
relatively constant pressure. Top and side views are shown. 
The side view is looking axially toward the carrier. The 
three dotted lines show single paths by which bubbles may rise 
through the sample volume. The paths are labeled by the 


maximum carrier signal pressure encountered. The exciter 
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TABLE 12 SAMPLE VOLUME SIZE CHARACTERISTICS FOR @ = 55° 


CARRIER 
PRESSURE a Dae AA yy AH 
‘(aB re MAX) (mm) = (mm) = (mm*), (mm) (mm) 
0 to -1 0.2 O20. Oee0 0.4 0.4 
-1 to -3 O52 0.44 0.24 0.6 One 
-3 to -6 0.5 Pe 1). 7 1.0 0.4 
-6 to -14 0.8 Biol Secon on 1.6 0.6 
-14 to -20 1.0 4.88 1.75 26 0.4 
-20 to -35 Ves 10.99 6.11 30 na 


uo 


ouNTOA eTdues sy} JO meTA PeTJTtTduts gst SInBTy 
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field is assumed constant and maximum over this small area. 
Table 13 depicts the time history of 100 um, 204m, and 9 um 
bubbles as they rise along these three paths. fe 1S "“elear 
that a 30 kHz bubble stays in the sample volume no longer than 
0.17 seconds. Though this time is enough for the bubble to 
be excited it is very short for the receiving equipment to 
detect the signal. If rise time was the only consideration, 
then any bubble smaller than 30 wm radius would be in the 
sample volume at least 1 second. Unfortunately, radiation 
force of the exciter and carrier push bubbles out of the 
Sample volume. The smaller the bubble, the more it 1s pushed 
and less likely to stay in the sample volume. To reduce this 
effect, the maximum pressure must be reduced. This, in turn, 


reduces the sample volume further. A compromise must be made. 


D. ELECTRONIC EQUIPMENT 

This section will detail the primary experimental setup 
in two major sections, drive signal generation and receiving 
equipment setup. The total system diagram is shown in Figure 
io, 

1. Drive Signal Generation 

The carrier signal was a 2.65 MHz sine wave. Because 

the sum and difference signals flank the carrier on a spectral 
plot, any anomalies near the carrier frequency could not be 


tolerated. This criterion led to the selection of the HP 
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TABLE 13. TRAVEL TIME AND CARRIER PRESSURE FOR BUBBLES 
TRAVELING ALONG THE PATHS SHOWN IN FIGURE 18 
TIME (mSEC) PRESSURE (dB re MAX) 

BUBBLE 100 wm 20 um 9 um PATH OGB -6 dB -20 GB 

0 0 0 -35t ~~ a= 
5 89 434 -35t -- 
22 398 1939 -35t 
28 521 2538 -20t 

oil 660 3216 -20t 

40 729 3552 -~14t 

50 912 4444 -14t 

57 1042 5077 - 6t 

68 1250 6091 - 3t 

74 1354 6598 - lt 

85 1563 7616 Ox — 6% -20% 
96 1771 8629 - 1) 

102 1875 9136 - 31 

Pe: 2083 10150 - 6! 

120 2213) 1o7se -14) 

130 2396 11675 -14) 

134 2465 12011 -20) 

142 2604 12688 -20} 

148 2727 oe ee ~354 

165 3036 14793 -~351 

170 3125 15288 -351 


t INCREASING PRESSURE, * MAXIMUM PRESSURE 
4 DECREASING PRESSURE 
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8640B Signal Generator which has very good spectral purity. 
C5575 could be driven to provide a maximum focal pressure of 
about one and a half atmospheres without using an amplifier. 
The voltage was measured with a HP 400E analog AC voltmeter 
or the oscilloscope. 

The excitation signal was originally generated by a 
HP 3314A Function Generator in the sweeping mode. This method 
produced bubble detection in most cases, but proved to be 
somewhat ambiguous in bubble size estimation. With the 
sweeping function generator, any given bubble's resonance 
frequency is broadcast once in the sweep interval. However, 
by broadcasting band limited white noise every frequency is 
present (though with randomly fluctuating amplitude), 
increasing the probability of detection. A HP 8904A 
Multifunction Synthesizer was used to generate the band 
limited white noise. An Amplifier Research model 50A15 
amplifier with a maximum power of 50 watts and a frequency 
range of 20 kHz to 15 MHz was used to amplify the noise. A 
passive low pass RC filter was used to eliminate frequency 
components above 500 _ kHz. This signal was monitored 
periodically on the spectrum analyzer to verify proper 
frequency coverage. An oscilloscope was used to verify the 
overall noisy quality of the signal. The exciter was driven 
at high power to ensure enough energy was present at each 
frequency to excite any bubble in the size range of interest 


into large amplitude resonant oscillations. 
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2. Receiving Equipment 

The receiving system consists of a common branch 
supplying three independent detection devices. The receiving 
transducer converts acoustic signals to electrical signals. 
This signal passes through a DC blocking capacitor and an 
amplifier with a constant gain of 27.7 dB over the range of 
interest. Three paths diverge from this point. The first 
goes directly to an oscilloscope. The second into the 1 
megaohm input to the HP 3585A Spectrum Analyzer. The third 
path leads to a HP 3561 Dynamic Signal Analyzer via a 
frequency mixer. This third path was often disconnected due 
to certain limitations of the mixer. The functions of these 
detection devices are discussed in this section and the 
results chapter. 

The HP 3585A Spectrum Analyzer is the heart of the 
development version of this bubble detector. The oscilloscope 
was useful, but not capable of yielding quantitative 
information due to the short duration of the signal and the 
overwhelming magnitude of the scattered carrier compared to 
the sum frequency signal. A frequency mixer was employed to 
take the high frequency sum and difference signals down to 
base band. Then the HP 3561 Dynamic Signal Analyzer could 
search the first 100 kHz to detect the larger, faster bubbles. 
This setup was both useful and frustrating. Neither of these 
instruments are continuously attentive to all frequencies. 


Bubbles were missed while these machines were busy with 
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calculations or searching other frequencies. Methods were 
hypothesized to correct this, but not implemented. These 
methods will be discussed with the conclusions. 

The oscilloscope would display the real time received 
Signal. Without bubbles in the sample volume, this display 
was a sine wave at the carrier frequency. During a bubble 
event, this display would increase dramatically in amplitude 
as the bubble would scatter the exciter and carrier, as well 
as, the sum and difference frequencies toward the receiver. 
A sharp increase in the oscilloscope trace amplitude would be 
interpreted as a scattering event. It was not possible to 
tell what type of scatterer caused the event due to the speed 
of the fluctuation. This was used as an indication that a 
bubble actually passed through the sample volume. 

The HP 3585A uses a frequency sweeping receiver. The 
receiver is sensitive to a narrow bandwidth that sweeps across 
the frequency range specified. The fastest sweep time is 0.2 
seconds. A constant signal will be recorded each sweep. An 
intermittent signal is recorded only if the sweep is covering 
that frequency coincidentally. By adjusting the frequency 
range to sweep both the sum and difference frequency, each 
bubble has two chances to be detected each sweep. Hence the 
frequency range for bubble detection was 2.25 MHz to 3.05 MHz. 
This permitted detecting bubbles of up to 400 kHz resonance 


frequency. 
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The HP 3561 uses a fast Fourier transform (FFT) 
technique. It records an input signal of at least 4 
milliseconds, and then computes and displays the frequency 
spectrum of that signal via a FFT. It does not record data 
during the computation. An alternative to this method is to 
take 40 consecutive records of 4 milliseconds, then analyze 
and display the whole batch. This method is both better and 
worse. A greater period of time is covered during the 
consecutive records, but a greater period is left uncovered 
during the computations. 

The mixer required for the use of the HP 3561 has two 
distinct disadvantages. First, it distorts the received 
Signal mildly in the region of the sum and difference 
frequencies as seen on the HP 3585A. Second, its own power 
limitations limit the power that can be simultaneously 
supplied to the carrier. Of these two, the distortion is the 


more serious problem. 


E. MODE OF OPERATION 

Many experiments were performed during the development of 
this system. Two major tests comprised the validation of the 
method to detect and size bubbles. The first involved timing 
a Single bubble as it rose through a known distance. The 
resonance frequency corresponding to the rise time was then 
compared with the resonance frequency corresponding to sum and 


difference frequencies indicated on the HP 3585A. The purpose 


oo 


of this test was to show that the sum or difference 
frequencies were being detected correctly, and that the sum 
or difference frequency gave an accurate and consistent 
indication of the bubble's resonance frequency. 

The other major validation test involved sending a single 
bubble through the sample volume and comparing the outputs of 
the HP 3585A, HP 3561 and oscilloscope. This test was 
designed to show that, at least for bubbles having resonance 
frequencies less than 100 kHz, the spectrum analyzer gave 
results consistent with those obtained from the HP 3561 Signal 
Analyzer. The two analyzers process signals differently. The 
Spectrum analyzer samples a given frequency component only 
once per sweep. The signal analyzer digitizes a complete 
window of data and then performs a Fourier transform on that 
data. In that sense, the signal analyzer processes all 
frequency components simultaneously. 

The two validation tests had one significant limitation. 
Bubbles smaller than 30 um radius cannot be easily seen by the 
naked eye to be timed and the HP 3561 cannot detect bubbles 
with resonance frequencies greater than 100 kHz. Therefore, 
resonance frequency determination of bubbles smaller than 30 
um radius (which also have resonance frequencies greater than 
100 kHz) cannot be verified by either method. In order to 
partially validate the system at the higher resonance 
frequencies, coincidence was used. If a very small radius 


bubble was detected by the system after the bubble generator 
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was momentarily energized, but no bubble was seen, then the 
system would be believed. This is not validation, but it is 
reasonable extrapolation. 

Other tests were performed with multiple bubbles in order 
to determine the systems response in this situation. Bubbles 
were readily detected, but the system display saturated 
quaickly . This effect will be explained in the results 


section. 
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IV. RESULTS 


The results of the tests discussed in the previous section 
are presented in this chapter. However, a brief discussion 
of the output display of the spectrum analyzer is presented 
firse The results of the single bubble detection tests, 
including data from the rise time resonance frequency 
verifications, follow. Tests of the dependence of the output 
signal level on the exciter and carrier pressure levels are 
then considered. Finally, the results of injecting solid 


particles into the sound field are discussed. 


A. TYPICAL OUTPUT DISPLAY 

The HP 3585A is the primary detection device for this 
system. This spectrum analyzer sweeps its frequency window 
through the range specified by the operator and displays a 
graph of amplitude versus frequency. Since the signal from 
a passing bubble is a fast acting transient phenomena, the 
most useful display mode for the purposes of this thesis is 
the "maximum hold" mode. In this mode the current display 
reflects the highest amplitude that has been present in a 
given frequency bin since the start of data acquisition. The 
use of this mode not only results in the display of the 


maximum signal detected but also the maximum noise. 
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In order to recognize a bubble detection, the display must 
be understood. Without bubbles, a small portion of the 
exciter and carrier signals are detected by the receiver. 
Figure 20 is a plot of the HP 3585A display of the received 
Signal under this case. Figure 21 is the same display after 
the mixer has been added to the receiver circuitry. Notice 
that the mixer has both increased the noise in the frequency 
range 2.25 MHz to 3.05 MHz and reduced the 2.65 MHz carrier 
Signal. Notice also the clean spectral nature of the carrier 
and the variations in the exciter level (0 ~ 400 kHz). The 
high amplitude signal from 4 to 5 MHz in Figure 21 is higher 
order modulation of the carrier signal with the received 
Signal. This noise is excusable because it is not in the 
bubble detection frequency range. 

In order to detect rapidly rising bubbles, the fastest 
sweep over the frequency region of interest was needed. This 
resulted in the 800 kHz span centered at 2.65 MHz being 
covered in a sweep time of 0.2 seconds. Any bubble with 
resonance frequency, fr, in the range 30 to 400 kHz would have 
both a sum (2.65 MHz + fp) and difference (2.65 MHz - fp) in 
this frequency band. Figure 22 shows the display of this 
region as it appears without bubbles (and without the mixer). 
Figure 23 shows the same frequency span after copious bubbles 
were allowed to rise through the sample volume for five 
minutes. A comparison of the two figures shows that the 
carrier level is at least 27 dB higher in Figure 23, 
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indicating that it was scattered significantly towards the 
receiver by the bubbles. Sum and difference frequencies were 
recorded for many bubbles of various sizes with a signal to 
noise ratio approaching 15 dB. Figure 24 is a superposition 
of the same type of plots for a higher carrier pressure. The 
bottom curve is without bubbles and without the mixer. The 
top curve is with bubbles but not the mixer. The middle curve 
has both bubbles and the mixer. As shown by the top curve, 
sum and difference frequencies are detected with a 10 to 20 
adB signal to noise ratio. These results are strong evidence 
which support the bubble detection goal. 

Figures 23 and 24 also illustrate system saturation. When 
the spectrum analyzer is in the maximum hold mode, the 
Signature of a bubble, once detected, is maintained on the 
display. This signal is now noise to future detections and 
the system becomes less sensitive to bubbles of size similar 
to those it has already detected. This saturation effect is 
not a problem for these experiments since they involved only 
one bubble at a time. However, it indicates that this 
processing arrangement is not well suited for multiple 
detections. 

B. SINGLE BUBBLE DETECTION AND RESONANCE FREQUENCY 

DETERMINATION 

Figures 25, 26 and 27 show three variations of single 
bubble detection and resonance frequency determination. In 


each case, both the sum and difference frequencies were 
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detected. Notice that the signal to noise level is a little 
different in each case. This difference is most likely due 
to the wide variety of carrier pressure available in the 
sample volume. 

The frequency of highest amplitude, other than the 
carrier, corresponds to the bubble's resonance frequency. 
Using the rise time method to independently determine 
resonance frequency verified the accuracy of the frequency 
determination. A large number of such comparisons were made. 
Figure 28 shows the agreement between the two mentiods: 
Perfect agreement would correspond to all of the data points 
falling on the line having a slope equal to one. The two 
methods agree to within a few kilohertz in most cases. A 
major source of error lies in timing the tiny rising bubbles. 

While performing these single bubble experiments it was 
noted that not every bubble was detected. There were some 
cases where bubbles passed through the sample volume 
sufficiently to show a transient on the oscilloscope but did 
not result in sum or difference frequency peaks on the HP 
3585A display. Since the oscilloscope verified the presence 
of a received signal, the HP 3585A was apparently inattentive 
to the frequencies of interest. The absence of these 
frequency peaks could have been due to: (1) the HP 3585A 
performing an internal calibration procedure that prohibited 
sampling while the bubble was in the sample volume, (2) the 


HP 3585A sweeping the wrong frequencies during a fast bubble 
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transient, (3) the bubble passing through a portion of the 
sample volume having too low a signal to noise level to permit 
detection, and (4) too much energy being scattered into the 
receiver causing the HP 3585A to overload. Some credence is 
given to all of these possiblities. Some displays would 
detect the sum or difference frequency but not both. Only 
alternative number (2) above makes sense in this case, since 
both signals should be present in approximately the same 
magnitude at the same time. 

The HP 3561 uses an FFT technique to analyze the spectral 
nature of the signal. However, the use of this instrument did 
not improve the probability of detection for several reasons. 
First, there is significant time over which the instrument is 
not sampling during computations, calibration and display. 
Second, the instrument's frequency range is limited to 0 to 
moom KHZ. Hence, it could never detect a 100 to 400 kHz 
bubble. Finally, it requires use of the mixer which 
inherently lowers the signal to noise ratio. Some low power 
detections might be lost in the noise. 

However, the HP 3561 was useful. First, it verified that 
FFT techniques can be employed in this problem. Second, it 
provided a different method to verify the accuracy of the 
resonance frequency determination. Figure 29 shows the 
Gisplay of the HP 3561 for the case without bubbles. This 
shows the noise floor as a function of frequency. Since the 


mixer demodulates the signal from a 2.65 MHz carrier, the 
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frequency scale labeled O Hz to 100 kHz actually represents 
Peeoo CO 2.75 MHZ. 

Figure 30 shows the detection and resonance frequency 
determination of a single bubble by the HP 3561. Using the 
time buffer mode, the HP 3561 continuously recorded 40 data 
records prior to processing. Ten of these data records were 
analyzed and plotted here. The first three and last six of 
the ten records show no hint of a bubble. Record four, on the 
other hand, strongly indicates the presence of a 40.25 kHz 
bubble. This figure indicates the transient nature of the 
bubble signals. The HP 3585A simultaneously detected the same 
bubble. Figure 31 is a copy of the HP 3585A display for this 
run. It shows the difference frequency strongly at 2.61 MHz, 
which represents a resonance frequency of 40.00 kHz. The sum 
frequency did not show up due to the high speed of this buble 
transient. In other words, the signal was gone before the 20 
milliseconds needed for the sweep to advance to the sum 
frequency had elapsed. The HP 3561 bears this out since each 
record represents only 4 milliseconds of data. Although two 
peaks of approximately the same amplitude are shown in Figure 
30 the one at 40.25 kHz is actually much larger. It appears 
to be the same amplitude because, it exceeded the maximum 
scale and was truncated. The level -87.13 dBV is, however, 


accurate. 
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C. VERIFICATION OF DUAL FREQUENCY METHOD PRESSURE 
DEPENDENCIES 


In order to confirm that the detected signals are in fact 
produced by the mechanism described by the dual frequency 
method, the peak pressure amplitudes of the detected signals 
were compared to the predictions of equation 2. The 
calculations were performed for many different bubble sizes 
using the transducer calibrations to determine the maximum 
Carrier pressure, the exciter pressure at the bubble resonance 
frequency and the pressure amplitude of the sum or difference 
frequency. Figure 32 is a comparison of the actual received 
level and maximum calculated sum or difference pressure at the 
receiver face divided by the product of carrier and exciter 
pressure. The expected values were given earlier in Table 4. 
None of the actual values exceeded the maximum expected 
values. The lower actual values are attributable to the fact 
that most bubbles encountered less than maximum carrier 
pressure. 

According to the theory of dual frequency bubble 
detection, the amplitude of the return signal depends on the 
pressure amplitudes of the carrier and excitation at the 
location of the bubble. In order to get an idea of what drive 
pressures were required, the exciter pressure was varied over 
a range of 7.36 dB with a constant carrier pressure amplitude 
of 102 kPa. Since the carrier pressure varies dramatically 
over short distances and the excitation pressure varies wildly 


with frequency, it would take an enormous amount of data to 
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Figure 32 Comparison of Actual Pressure Relationship 
with Maximum Expected 
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verify equation 2 with this implementation. In spite of this, 
Figure 33 shows that the signal level varies with the exciter 
pressure. The exciter pressure is lowest for frequencies of 
about 140 kHz. Therefore, 140 kHz bubbles are detected with 
the weakest signal. The main reason for the curves not being 
flat is the variations of excitation pressure with frequency. 
Another factor is that smaller bubbles naturally have smaller 
Signals for the same excitation. 

Figure 34 shows a similar family of curves. This time the 
exciter was held constant at the maximum value. The carrier 
pressure was varied by a factor of 2 (6dB). Again, the bubble 
Signal decreased as carrier pressure decreased. It is clear 
that for better signal to noise ratios, the carrier and 


exciter pressures should be as high as possible. 


D. DISCRIMINATION OF BUBBLES FROM SOLID PARTICLES 

Dry grains of sand, ranging in size from 1 um to 1 mm 
diameter were dropped into the water over the sample region. 
The grains rarely entered the sample volume, but when one did 
it almost always produced a signal similar to that produced 
by bubbles. These signals were attributed to the presence of 
entrained gas pockets, some of which were visible, on the dry 
sand. It was reasoned that wet sand would not produce the 
same type of signals. Therefore, pre-wetted sand was more 
carefully injected into the sample volume. Only occasionally 


would an extraneous bubble indication be noted under these 
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Figure 33 Received Signal Variations with Changing 
EX@heer ERessuire 
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Received Voltage Level (dBv) 
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Figure 34 Received Signal Variation with Changing 
Carrier Pressure 
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circumstances. This was thought to be due to an infrequent 
gas bubble in the slurry. The most common response Cfo tne 
system to a solid is an increase in the received carrier and 
exciter signals. As long as the pressure of these signals is 
less than that required for medium non-linearities, such as 
the parametric effect, no sum or difference frequencies should 
be detected. 

Hollow glass beads were also injected into the sample 
volume. The response was identical to that for gas bubbles. 
This response semed reasonable since the glass beads were at 
least 90 percent gas by volume. They floated rapidly to the 
surface. The results of these measurements indicate that the 
dual frequency method is capable of detecting bubbles as well 
as other partially gaseous matter. However, it is insensitive 


to non-gaseous scatterers. 


E. DISCUSSION 

ise to this point, many facts and figures have been 
presented, but the full implications have not been discussed. 
This section will review the evidence presented for this 
implementation of the dual frequency method for detecting and 
sizing bubbles. 

The first requirement was to detect bubbles. This goal 
clearly was accomplished. Bubbles were detected through the 
range of interest, though the majority of the events were less 


imam 100 kHz. Data to support this was taken with single 
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bubbles and lots of bubbles, free rising bubbles and convected 
bubbles, varied excitation and carrier pressure levels. By 
Going this, not only were bubbles detected, but some 
limitations of bubble detectability were discovered. These 
will be highlighted later in this section. 

The ability to directly determine the resonance frequency 
of a bubble has been shown. Two independent signal processors 
were used to show that this technique was not based on a 
peculiarity of one type of processor. The resonance frequency 
determination follows the simple rule that the bubble has its 
maximum response at resonance. Therefore, the maximum 
received signal in the detection band corresponds to the 
resonance frequency. The rise time method of resonance 
frequency estimation supported the dual frequency method 
results. The agreement in the results of the two processors 
and the rise time method was best when detected with high 
Signal to noise. This was not surprising since low signal to 
noise was often due to poor bubble signal timing with respect 
to the frequency sweep for the HP 3585A or the signal data 
record for the HP 3561. Strong signal to noise indicated that 
the bubble was properly detected while in the processor's 
analysis window. In these cases, this dual frequency method 
determines the resonance frequency to within the signal 
processor's analysis error. 

The tests with solid particles indicated that this method 


can determine the resonance frequency of solid borne or 


1) 


membrane enclosed gas volumes. The accuracy of these 
measurements has not been shown, but the presence of a signal 
very similar to that of a pure bubble for these gas volumes 
with solids suggests the same ability to correctly determine 
resonance frequencies. A solid particle without gas was shown 
to lack the sum and difference frequency signals. Therefore 
non-gaseous scatterers are easily distinguishable from various 
types of gaseous scatterers having resonance frequencies in 
the range of interest. This natural ability to classify a 
bubble by its acoustic resonance parameters has far reaching 
implications for the study of bubbles in both clean and dirty 
environments. 

During the bubble experiments, a wide variety of signal 
to noise levels were noted. Some bubbles transited the sample 
volume without any indication of its presence on the 
processor's display. Others were detected with more than 20 
@B of signal to noise. Much of this was probabilistic since 
the bubble's path through the sample volume and signal timing 
to the signal processor were not controllable. This random 
positioning and timing meant that any single run with any 
Single bubble could not be directly related to other runs. 
To compare one test setup with another many data runs were 
taken. The maximum signal to noise cases were considered to 
have had the most optimal sample volume positioning and 
timing. Thus, the differences between these runs would be due 


to the actual parameter changes between test setups. An 


ne a 


example of a test setup change is the exciter pressure 
variations previously discussed and shown in Figure 33. To 
verify the relationship of bubble signal pressure to exciter 
pressure many bubble runs were made at each pressure. The 
maximum bubble signal at each pressure was considered to be 
optimum in position and timing. 

The sum frequency pressure relationship given by equation 
2 was tested by doing these comparisons. To establish a 
threshold of detection the whole system must be understood. 
For this system, the detection band noise level was about -100 
adBv. To obtain 10 @B signal to noise, a signal of -90 dBv is 
needed at the instrument. Subtracting the 27.7 dB amplifier 
gain yields a necessary -117.7 daBv out of the receiver. It 
is easy to calculate the pressure required at the receiver 
face to generate this voltage. Using equation 2, the required 
excitation pressure can be found for each frequency given the 
Carrier pressure. Table 14 shows the results of these 
calculations. 

The pressure produced by the exciter was well above the 
calculated values of Table 14. The carrier pressure was also 
higher than the value assumed for this table. These higher 
pressures produced the data with 20 dB signal to noise. The 
weaker Signals were due to the losses induced by poor 
positioning and timing. The numbers in Table 14 are valid for 


any system since they are simply pressure relationships. The 
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TABLE 14 
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EXCITER PRESSURE REQUIRED FOR 10 dB SIGNAL TO NOISE 
GIVEN r = 


30K 10° Pa 
Exciter 
Pressure FREQ. 
(Pa (peak) ) (KHz) 
29.6 50 
68.0 100 
BEES 150 
184 200 
254 250 
274 300 
319 350 
374 400 


receiver output voltage would, of course, vary with different 
receivers. This makes this analysis specific to this system. 

If a system could be designed to produce these exciter and 
carrier pressures and receiver coverage over a broad sample 
volume, it would detect all bubbles with about 10 dB signal 
to noise ratio. This assumes that a larger sample volume 
would alleviate the critical timing problem since the bubble 


would remain in the sample volume longer. 
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V. CONCLUSIONS 


The goals of this bubble detection system have been 
accomplished. The system has demonstrated the ability to 
detect bubbles with resonance frequencies from 30 to 400 kHz 
when they entered the sample volume. The resonance 
frequencies of these bubbles were accurately determined by the 
system. The ability to distinguish between bubbles and non- 
gaseous particles was shown. Parameters for creating a system 
with a high probability of bubble detection were accumulated. 
The ability to determine the resonance frequency of scatterers 
containing a small gas volume was also noted. 

With the achievement of these goals, the dual frequency 
method for bubble detection and resonance frequency 
determination becomes a powerful tool. It can be used to 
analyze bubble populations or to study the dynamics of single 
bubbles under the influence of stress. Studies may be done 
using this technique to study the effect of varying a fluid's 
Viscosity or surface tension on a bubbles damping constant. 
There are a number of areas for application of this technique. 
The next section recommends items to consider to improve the 


performance of a dual frequency method bubble detector. 
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A. RECOMMENDED IMPROVEMENTS FOR A CONTINUOUS BUBBLE DETECTION 


SYSTEM 


Many items have been mentioned that effect the chances of 


detecting a bubble with this system. Some improvements must 


be made to achieve an ideal system. Here is a summary of 


areas of concern. 


* Sample Volume 


* Bubble 


x Receiving Equipment 


x Environment 


Concerns 


too small 
too variable 
too weak 


path and time in sample volume 


needs to sample continuously 
needs large dynamic range (to 
detect small sum and 
difference frequency signals 
in the presence of large 
scattered carrier and exciter 
signals) 

data output management 


noise level in detection band 
number of bubbles: saturation 
effect and signal extinction 
vibrations: disturb sample 
volume alignment 

flow rate: speed of bubble 
through sample volume 
quantity of suspended solid 
scatterers 


1. Sample Volume Related Improvements 


A bubble must pass through the sample volume to be 


detected. A large sample volume permits easier coverage of 


reasonable volumes. The smaller the sample volume is with 


respect to the total systen, 


1G 


the lower the probability of 


detection. The implementation tested here had virtually the 
smallest possible sample volume. 

Variations in the carrier and excitation pressures and 
receiver efficiency within the sample volume allow for highly 
variable received signals when a bubble is present. If the 
lowest possible received signal is less than the noise, 
detection opportunities will be missed. Driving the exciter 
with band limited noise in the frequency range of interest 
provides broad frequency coverage, with only minor temporal 
variations. The drive current and exciter source strength are 
functions of frequency. The drive current should be 
manipulated to excite bubbles equitably, to give each size 
bubble the energy to obtain the same signal to noise. 

The receiver's efficiency may vary with frequency and 
bubble location. These variations are strictly functions of 
the transducer. The transducer should be designed or selected 
to minimize these variations in the frequency range and 
spatial region of interest. The carrier transducer is 
monofrequency. Hence, only stability of output power and beam 
pattern are of concern. The overall goal is to obtain a good 
Sized volume of constant and adequate, but not excessive, 
acoustic pressure. If this is accomplished, many of the 


problems associated with a small sample volume disappear. 


na key 


2. Bubble Related Concerns 

The path of travel varies from bubble to bubble, due 
to a variety of forces acting on the bubble. Gravity, 
buoyancy, drag and radiation forces are the key players. Of 
these, radiation force may be the worst. This acoustic 
phenomena pushes the bubble when the acoustic pressure is 
high. The resulting motion drives the bubble toward regions 
of lower carrier and excitation pressure. Keeping the 
pressure just below the level that causes significant motion 
is important to solving this problem. Having a relatively 
large sample volume would help regardless of the cause of the 
bubbles motion. 

The amount of time a bubble spends in the sample 
volume impacts its detectability. It is a function of the 
sample volume, bubble speed and path. These items need to be 
manipulated to achieve an ideal length of time in the sample 
volume. Too little time results in missed detections. Too 
much time will limit the rate at which the system can 
distinguish between individual bubbles. 

3. Receiving Equipment Improvements 

The receiver must be alert to detect an event. The 
equipment assembled for this experiment was useful, but not 
ideal for the job. The received signal needs to be 
continuously monitored over the frequency range of interest. 
Only the sum or difference band is actually needed, but both 


may be wanted for verification. The data needs to be 
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presented in a format in which new detections are not effected 
by previous ones. A continuous hard copy spectrogram is a 
good way to handle large amounts of data. 

During a bubble detection, many frequency signals are 
received. Unfortunately, those of the most interest are of 
low level. Large dynamic range is needed in the receiver. 
This allows it to measure the low level signals without 
overloading due to the high level scattered signals. 
Filtering the return signal can help this problem. Band pass 
filtering to accept only the detection band followed by notch 
filtering to eliminate the carrier would eliminate the 
unwanted sSignals' overbearance. 

Correcting the overall alertness and relative dynamic 
range of the receiving equipment would be a major improvement 
in this system. This would require a system designed for this 
purpose. Most off the shelf equipment does not have the 
frequency range or data volume capability needed for a true 
bubble detection system. 

4. Environment Related Concerns 

The environment of the test apparatus has a large 
impact on its effectiveness. Vibrations change the bubbles 
path, as well as, the sample volume alignment. The quantity 
of scatterers throttles the amount of power available to 
excite a bubble and effects the propagation of the bubble's 
Signal to the receiver. The ambient noise in the detection 


band is a factor in detectability. The environment was not 


dO 


a major problem in this laboratory experiment. It has the 
possibility of causing problems, but they are minor compared 


to the previous three categories. 
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APPENDIX 


EASY REFERENCE TABLE 
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BUBBLE RES. FREQ DAMPING _p 1 TIME TO 
RADIUS, a (um) _ (kHz) CONSTANT _p;pp_ Pa RISE 1" (SEC) 
100 Bal, 4 05075 me.3 x 10° iL ge 
70 44.6 jl Osseoe x 107 2.55 
50 Dod 0.095" le 74 x TO A261 
40 ee Coins 4) x 10i- 6.95 
30 102.6 Geolismeece.79 x 10° 12-0 
20 152.9 Ors 5.25 x 10° 26.50 
10 303.8 0.16 mG x Lo 105 
9 32929 0.16 oA ye TO” 129 
8 380.4 0.16 2 sk 0” 163 
7 435.8 Ge 16 iL git 9% alae Pale 
6 510.8 0 fale 28 x 10” -- 
5 618.0 Neale 1.05 x 10°” 416 
EXCITER CARRIER RECEIVER 
FREQ PpEAK Vo Peeak FREQ MG Pa) 
(KHz ) (Pa) (Pa) (MHZ) 

50 115 0.052 FO 2.25 0.0000982 
100 BEM 0.165 2242 Page: 0.0000969 
150 244 0.41 5570.0 £2.35 0.0000955 
200 370 0.55 FAID40°. 2-4 0.0000951 
250 515 0.94 12800. 2.45 0.0000939 
300 980 AL. Ie, 14940. 2S 0.0000936 
350 1000 1.40 18250. 2.55 0.0000932 
400 1190 1.98 DiESVAON 2G 0.0000934 

2.08 28260. 2.65 0.0000932 
2.55 SAGA Or DG 0.0000937 
3.90 51800. DATS 0.0000948 
4.0 47080. 288 0.0000947 
4.05 55020. 2.85 0.0000967 
5.0 66580. Doe 0.0000966 
5.2 70650. 2.95 0.0000982 
6.0 81540. 3.0 0.0000974 
7.0 88080. 3.05 0.0000972 
7.67 102400. ela 0.0000966 
8.0 105300. 
8.2 111400. 
9.0 115300. 
ORO 133200). 
10.35 140600. 
TI 148900. 
T2029 175300. 
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